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ABSTRACT 


The  Lineate  Imaging  Near-ultraviolet  Speetrometer  (LINUS)  uses  three  eomplex 
lens  systems  to  focus  an  image  from  distances  on  the  order  of  several  kilometers  onto  the 
image  intensifier  of  an  ultraviolet  camera.  These  images  can  then  be  analyzed  to 
characterize  the  atmospheric  distribution  and  concentration  of  sulfur  dioxide  (SO2).  The 
lenses  purchased  for  LINUS  were  corrected  for  spherical  aberrations  but  due  to  the  lack 
of  detailed  knowledge  about  the  lenses  their  chromatic  aberrations  could  not  be  readily 
predicted.  The  project  presented  in  this  thesis  was  performed  with  the  goal  of 
experimentally  quantifying  the  chromatic  aberrations  of  each  of  LINUS’s  lens  systems. 

The  matrix  method  of  representing  paraxial  optical  systems  was  used  to  determine 
relationships  between  object  and  image  distances  for  different  lens  systems.  These 
relationships  were  manipulated  to  give  equations  for  the  matrix  elements  of  the  lens 
system  in  terms  of  readily  measurable  parameters.  Once  the  matrix  elements  are  known, 
all  of  the  cardinal  points  can  be  readily  predicted.  This  method  will,  in  theory,  quantify 
the  chromatic  aberrations  of  each  lens  system.  The  method  was  validated  with 
simulations  and  measurements  taken  on  a  lens  of  known  focal  length.  Finally,  the 
LINUS  lens  systems  were  characterized  at  220,  300,  334,  and  370  nm. 
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I.  INTRODUCTION 


A.  MOTIVATION  AND  OBJECTIVES 

Ultraviolet  (UV)  and  near-ultraviolet  speetroseopy  researeh  has  been  ongoing  at 
the  Naval  Postgraduate  Sehool  sinee  1990.  The  first  sueh  instrument  developed  here  was 
the  Middle  Ultraviolet  SpecTrograph  for  Analysis  of  Nitrogen  Gases  (MUSTANG),  a 
roeket-borne  instrument  intended  to  analyze  emission  line  speetra  in  the  earth’s 
ionosphere.  MUSTANG  was  later  upgraded  from  a  eonventional  speetrograph  to  an 
imaging  speetrometer  resulting  in  the  Dual  Use  Ultraviolet  Imaging  Speetrometer 
(DUUVIS).  DUUVIS  was  tested  early  in  1997.  However,  the  images  that  it  produeed 
were  not  of  aeeeptable  quality,  and  a  new  instrument  was  needed.  The  next  step  toward 
hyperspeetral  imaging  was  taken  with  an  entirely  new  instrument,  the  NPS  Ultraviolet 
Imaging  Speetrometer  (NUVIS).  NUVIS  beeame  operational  in  the  fall  of  1997  and  was 
field-tested  using  ehemical  effluent  plumes  from  industrial  smokestaeks.  Results  from 
this  testing,  analyzed  by  LT  Stephen  Marino,  demonstrated  the  ability  of  the  instrument 
to  eharaeterize  sulfur  dioxide  (SO2)  distribution  and  eoneentration  in  industrial  ehemical 
plumes.  LT  Marino’s  work  also  revealed  several  limitations  inherent  in  NUVIS’s  design 
[1].  One  of  the  more  severe  limitations  of  NUVIS  was  its  inability  to  observe  any  spectra 
other  than  that  of  SO2.  Rather  than  attempt  to  upgrade  NUVIS,  it  was  decided  to  build  a 
new  instrument  capable  of  hyperspeetral  imaging  across  the  entire  200-400nm  ultraviolet 
spectral  region.  Thus,  the  design  for  LINUS,  the  Lineate  Imaging  Near-Ultraviolet 
Spectrometer  was  conceived.  [1] 

Traditionally,  the  unavailability  of  acceptable  refractive  materials  with  low  UV 
absorption  properties  has  required  that  UV  instruments  consist  of  reflective  optics. 
However,  LINUS’  design  is  based  on  both  reflective  and  refractive  optics  due  to  the 
recent  development  of  low-absorption  UV  glass.  This  has  permitted  the  design  of  an 
aplanatic  instrument,  one  whose  optics  are  corrected  for  several  basic  aberrations. 

The  research  described  in  this  thesis  has  focused  on  the  characterization  of  the 
complex  lens  systems  used  in  the  LINUS  instrument.  Each  lens  system  was 
characterized  at  several  different  wavelengths  to  ensure  accurate  operation  at 


1 


wavelengths  other  than  that  eharacteristic  of  SO2.  Several  steps  were  taken  in  the 
solution  of  this  particular  problem.  First,  the  method  for  determining  the  lens  matrix  for 
each  lens  system  was  developed.  Next,  a  series  of ‘thought  experiments’  were  conducted 
to  verify  the  utility  of  the  selected  method.  Once  the  method  had  been  verified  it  was 
tested  using  a  thin  lens  and  visible  light  in  the  laboratory.  Finally,  data  was  taken  and 
analyzed  using  the  actual  LINUS  lens  systems. 


B,  OUTLINE 

This  thesis  is  comprised  of  five  chapters  and  two  appendices.  Chapter  II  contains 
a  brief  description  of  the  ray  transfer  matrix  method  of  characterizing  optical  systems  and 
describes  the  method  used  to  characterize  simple  lens  systems.  It  also  contains  the  results 
of  the  ‘thought  experiments’  used  to  validate  the  method.  Chapter  III  outlines  and 
verifies  the  method,  which  can  be  used  to  analyze  the  characteristics  of  complex  lens 
systems  in  a  manner  very  similar  to  that  used  for  simple  lens  systems  in  Chapter  II. 
Chapter  IV  presents  the  measurement,  analysis,  and  characterization  of  each  LINUS  lens 
system  at  selected  wavelengths  across  the  near-UV  spectral  region.  Conclusions  and 
recommendations  for  future  work  are  discussed  in  Chapter  V.  The  appendices  contain 
both  raw  experimental  data,  and  various  Matlab  programs  written  to  perform  the  data 
analysis. 
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II.  DETERMINATION  OF  MATRIX  ELEMENTS  FOR  SIMPLE 

OPTICAL  SYSTEMS 


A,  THE  MATRIX  METHOD 

Approximately  25  years  ago  several  innovative  researehers  developed  a  simple 
matrix-algebraie  approaeh  for  the  geometrieal  optieal  eharaeterization  of  paraxial 
systems. [2]  Due  to  its  ease  of  use,  this  method  was  the  one  ehosen  as  the  basis  for  this 
thesis  researeh.  The  basie  eoneepts  of  this  method  are  presented  here  as  an  aid  to 
understanding  how  this  method  was  employed  in  the  determination  of  thin  and  eomplex 
lens  system  eharaeteristies. 

1,  Paraxial  Optics 

Use  of  the  matrix  method  hinges  on  the  assumption  that  the  optieal  system  being 
analyzed  is  paraxial.  The  paraxial  behavior  regime  applies  whenever  the  angle  between 
the  system’s  optieal  axis  and  the  ray  of  interest  (a,  in  Fig.  1)  is  small.  This  allows  use  of 
the  small  angle  approximations  (sin(a)  ~  a,  tan(a)  ~  a,  and  oos(a)  ~  1)  when  traeing  the 
path  of  the  ray  through  the  optieal  system.  A  simple  paraxial  system  is  shown  in  Figure 
1. 


Thick  Lens 


2.  Matrix  Representation  of  Geometrical  Optics 

Each  of  the  three  arrows  seen  above  represents  a  translation  of  the  light  from  one 
point  in  space  to  another.  In  each  of  the  three  regions,  the  ray  follows  a  straight  line- 
segment  that  makes  a  constant  angle  between  the  ray  and  the  optical  axis.  At  points 
where  the  ray  meets  a  lens  face  or  other  discontinuity  in  refractive  index,  Snell’s  law 
dictates  that  the  ray  change  its  angle  with  respect  to  the  optical  axis.  This  constitutes  a 
point  operation. 


3 


First  one  must  find  a  way  to  define  any  point  in  the  ray’s  trajectory.  The  easiest 
way  to  accomplish  this  task  is  to  use  an  angle  and  height  format.  Height  of  the  ray  (y)  is 
measured  from  the  optical  axis.  The  first  ray  in  Figure  1  has  an  origin,  the  input  plane, 
and  an  ending  point,  the  first  face  of  the  lens.  The  point  of  origin  is  written  (ai,  yi),  and 
the  ending  point  is,  similarly  (ai,  yi)-  Assume  that  this  ray  has  traveled  a  distance  L, 
along  the  optical  axis.  The  actual  distance  that  this  ray  has  translated  is  L-cos(a)  but 
since  a  is  so  small  the  paraxial  approximation  applies  and  the  distance  translated  is  very 
closely  approximated  by  L.  Similarly,  the  ray  has  moved  away  from  the  optical  axis  by 
the  amount,  L-tan(a)  which  can  be  approximated  as  L-a.  The  following  relationships  can 
be  determined  simply  by  inspection  of  the  ray  in  question. 

a2  =  ai,  and 

yi  =  yi  +  aiL. 

These  equations  can  be  written  in  matrix  form  as  follows: 


Before  the  ray  undergoes  its  second  translation  the  angle  that  it  makes  with  the 
optical  axis  changes,  or  refracts,  with  no  change  in  height.  In  this  case  y2  is  obviously 
equal  to  y3.  The  change  in  angle  is  governed  by  Snell’s  law,  which,  in  its  paraxial  form, 
gives  the  following  relationship. 

«3=  ^  •  — -1  •T2+  —  «2- 

\RJ  \n  J  yn  ) 


Where  n  and  n’  are  the  indices  of  refraction  on  the  left  and  right  sides  of  the  interface  and 
R  is  the  radius  of  curvature  of  the  interface.  When  combined  with  the  fact  that  height 
stays  constant  for  this  point  operation  the  equations  can  be  written  in  the  matrix  form. 


.[3,  pp.  66-67] 

aij 
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Any  paraxial  optical  system  may  then  be  modeled  as  a  sequenee  of  physieal 
proeesses,  sueh  as  translations  and  refraetions.  The  individual  matriees  are  multiplied 
together  in  the  eorreet  operational  order  to  prediet  ray  paths  through  the  entire  system. 
For  example,  the  lens  shown  in  Figure  1  may  be  eharaeterized  by  a  matrix  as  follows:  let 
the  first  and  seeond  refraetions  be  represented  by  matrices  Ri  and  R2,  and  the  translation 
through  the  body  of  the  lens  be  represented  by  a  matrix  T,  where  Ri  and  R2  are  of  the 
form  M2  and  T  is  of  the  form  Mi.  The  system  matrix  for  the  lens,  whieh  we  will  eall 
Miens  is  defined  as:  Miens  =  R2’T-Ri.  This  matrix  eontains  all  of  the  information 
neeessary  to  deseribe  the  lens.  The  simplest  lens  matrix  is  assoeiated  with  a  perfeet  thin 
lens  of  foeal  length/  is  as  follows: 


= 


.[3,p.  70] 

3.  The  System  Ray  Transfer  Matrix  and  Its  Elements 


1  0 
1 


- 1 

/ 


A  eomplete  optieal  system  is  deseribed  by  an  input  plane,  where  rays  originate,  a 
sequenee  of  optieal  operations,  and  an  output  plane,  where  the  ray  pattern  is  analyzed. 
Thus  the  eomplete  system  representation  contains  two  additional  translation  matrices  that 
represent  the  rays’  path  from  the  system  input  plane  to  the  first  interfaee,  and  from  the 
last  interfaee  to  the  system  output  plane.  If  translation  from  the  origin  to  the  first 
interfaee  is  represented  by  Tin  and  the  translation  from  the  seeond  interfaee  to  the  output 
plane  is  represented  by  Tout,  the  system  matrix  (Msys)  for  the  system  shown  in  Figure  1 
will  be:  Msys  -  Tout'  R2'T-RrTin.  The  optieal  system  matrix  elements  are  constrained 
values  for  eertain  applieations,  sueh  as  fmite-eonjugate  imaging,  and  ean  be  used  to 
determine  the  loeation  of  eardinal  points. 


A  general  system  matrix  will  be  of  the  form: 


M  = 


A  B 
C  D 


When  a  finite-eonjugate  image  is  being  formed  at  the  output  plane,  matrix 
element  B  =  0,  and  A  is  the  linear,  or  transverse,  magnifieation  of  the  system.[3,  pp.  72- 
73]  These  faets  have  been  used  to  develop  a  method  to  experimentally  determine  the 

5 


system  matrix  for  real  lenses.  Onee  the  matrix  elements  are  known  for  a  given  lens,  it  is 
relatively  simple  to  determine  the  lenses  behavioral  eharaeteristies  for  various 
applieations. 

B,  EXPERIMENTAL  DETERMINATION  OF  MATRIX  ELEMENTS 

Given  an  unknown  lens  (or  eomplex  lens  system),  a  general  system  matrix  ean  be 
eonstrueted  by  ehoosing  arbitrary  input  and  output  planes  and  performing  finite- 
eonjugate  imaging  experiments.  An  initial  ray  translation  from  the  input  plane,  or  objeet, 
to  the  lens  by  a  distanee  So  is  followed  by  transmission  through  the  lens  (represented  by 
the  lens  matrix).  Finally,  there  is  translation  to  the  image  plane  loeated  a  distanee  s,  from 
the  lens.  The  matrix  for  this  complete  system  will  be. 


1 

to 

1 

-\ 

a  (d 

a  +  ys^ 

C  D' 

0  1 

y  5 

0  1  j 

y 

d  +  ys^ 

A,  B,  C,  and  D  are  the  system  matrix  elements.  As  stated  above,  for  finite-conjugate 
imaging  B  =  0,  and  A  =  linear  magnification  (m).  a,  ji,  y,  and  5  are  the  lens  matrix 
elements,  from  which  the  cardinal  points  may  be  determined.  When  the  above  system  is 
solved  for  the  lens  matrix  and  the  appropriate  substitution  made  for  the  A  and  B  elements 
the  following  system  of  equations  is  obtained. 


a  fd 

[1 

-1 

~A 

B~ 

"1 

y  5 

0 

1 

C 

D 

'  0 

1 

[1 

m 

o' 

"1 

0 

1 

C 

D 

0 

1 

m  -  Cs-  -Ds-  -{m-  Cs^  )s^ 

C  D-Cs^ 

Each  lens  matrix  element  can  now  be  expressed  as  a  function  of  the  system  matrix 
elements,  image  distance  {Si),  object  distance  {So),  and  linear  magnification  (m): 

a  =  m-  Cs., 

/3  =  -Ds,-{m-Cs,)s^, 
y  =  C,  and 
d  =  D-Cs^. 
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The  next  step  in  the  development  of  an  experimental  procedure  involves  reducing 
the  number  of  variables  in  the  equations  above  to  allow  expression  of  the  lens  matrix 
elements  in  terms  of  readily  measurable  quantities.  Since  each  of  the  lens  matrix 
elements  is  necessarily  constant  their  derivatives  with  respect  to  So  and  Si  must  equal  zero. 
Image  distance  (si)  is  assumed  to  be  the  independent  variable  and  object  distance  (^o)  and 
linear  magnification  (m)  the  dependent  variables  for  the  purposes  of  these  calculations. 
Then, 


da 

dSi 


dm  „  dm  .  , 

- L  - - y  =  0,  so  that 

ds,  ds, 


r  = 


dm 

ds, 


and  a  =  m 


Next,  the  d  equation  is  solved  for  D  and  substituted  into  the  fd  equation,  giving 


f3  =  -{d  +  Ys^)s.  -{m- rSi)s^  =  -ds.  -  ms^, 

dB  _  ds„  dm 

— ^  =  -d  -m  — ^  -  s^ - =  0, 

dsi  ds^  ds. 


„  ds^  dm  , 

d  =  -m — -- s„ - ,  and 

ds,  ds, 


P  -m 


ds^ 

s  — - —  s 
‘  ds,  ° 


dm 

ds. 


Each  of  the  lens  matrix  elements  has  been  expressed  in  terms  of  variables  that  can 
be  measured.  In  order  to  verily  these  equations  a  ‘thought  experiment’  will  be  performed 
in  two  steps.  These  expressions  will  first  be  validated  algebraically  using  the  familiar 
Gaussian  thin  lens  formula.  After  confidence  has  been  gained  in  the  mathematics  behind 
this  method,  a  set  of  data  will  be  generated  using  the  thin  lens  formulae  and  evaluated  to 
test  the  methods  utility. 
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C.  THE  THOUGHT  EXPERIMENT 

1.  Algebraic  Validation  of  the  Procedure 

This  portion  of  the  thought  experiment  involved  solving  the  thin  lens  formula  for 
eaeh  of  the  variables  in  the  equations  derived  above  for  the  lens  matrix  elements.  After 
the  variables  were  determined,  they  were  substituted  into  the  matrix  element  equations  to 
verify  that  the  expected  thin  lens  matrix  resulted.  The  thin  lens  equations  are, 

111 

—  =  — I - ,  and 

f 

m  =  -^. 


Expressions  are  needed  for  Si,  So,  m,  and  the  derivatives  of  So  and  m  with  respect  to 
Si.  Before  the  magnification  derivative  can  be  taken  however,  m  should  be  cast  only  in 
terms  of  st: 


s,  = 


fio 


s„  =  = — ,  and 


So-f 

_  A- 

A 


A-/)/-A_  -/ 
A-/)^  A-/)^’ 


dm  _  fs^ if -2s.)-  ifs.  -  sf  )/  _  1 

ds,  {fs,f  /■ 

These  expressions  will  be  used  to  show  that  each  lens  matrix  element  equation 
gives  the  expected  result  for  a  thin  lens. 


dm  s,. 

a  =  m-  s, - =  — --  s, 

ds,  s„ 


f. 


A 


/ 


A 
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P  =  m 


dSo 

S-  — -  —  s 
'  ds. 


dm 

ds. 


/•2 


-  S„ 


/ 


rZ  I 

—S- 1  s 

iJ  o 


2  1  ^ 


2  2 


V 


-  s„ 


f 


ss 


s„s. 


=  ^  +  s,-^  =  s^  +  s,-(s^  +  sp^0, 
s  s.  t 

O  I  J 


dm  1  , 

Y  =  —  = - ,  and 

ds,  f 


_  ds^  dm  s, 
o  -  -m  — -  -  s„ =  — 


f 


ds,  °  ds,  s^[{s,-f)  ) 


-f 


2  A 


/ 


—  s  s  s  s(s.  —  f)s 

J  O  _|_  o  _  o  _  o  _  o\  I  J  /  _  o  _  1 

^2  2"^  “  “  “  ~ 

\f^. 


f  f  s,  fs. 


Therefore, 


1  o' 

a  P 

_ 

1  , 

y 

—  1 

L  /  J 

as  expeeted.  This  indieates  that  the  proeedure  is  valid. 

2.  Utility  Verification 

As  previously  stated  the  purpose  of  this  portion  of  the  thought  experiment  is  to 
determine  if  the  method  of  determining  the  lens  matrix  elements  deseribed  above  will  be 
praetieally  aehievable.  The  following  simulation  was  performed:  artifieial  data  were 
generated  in  Mierosoft  Exeel  using  the  ideal  thin  lens  formula  for  a  hypothetieal  5  em 
positive  thin  lens.  Image  distanee  was  ehosen  as  the  independent  variable  and  assigned 
values  from  10  to  30  em  in  1  em  inerements.  Object  distance  was  calculated  from  the 
ideal  focal  length  and  image  distance  values.  Linear  magnification  was  then  calculated 
for  each  image  and  object  distance  set.  The  derivatives  of  magnification  and  image 
distance  with  respect  to  object  distance  were  estimated  by  calculating  the  differences 
between  magnification,  image  distance,  and  object  distance  as  shown  below. 


dm  Am  , 

- « - ,  and 

dSo 


ds.  Ay. 
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These  ratios  were  then  plotted  against  an  average  value  for  object  distance.  Data 
were  imported  into  Matlab  and  least-squares  fitted  to  a  user-defined  function  using  the 
Tsqcurvefit’  function.  The  user-defined  functions  were  defined  by  taking  the  appropriate 
derivatives  of  the  thin  lens  equation  and  are  as  follows: 


- and 

f 


Each  function  is  fitted  to  the  data,  which  is  formatted  in  (x,y)  pairs,  by 
determining  the  best  value  for  the  focal  length  f.  After  the  focal  length  is  determined, 
values  for  image  distance,  magnification,  and  the  derivatives  of  object  distance  and 
magnification  with  respect  to  image  distance  are  calculated.  These  values  are  then  used 
to  calculate  the  lens  matrix  elements  using  the  relationships  derived  earlier.  The 
simulated  lens  matrix  is  then  constructed  from  their  mean  values.  To  ensure  that  the 
curve  generated  is  a  reasonable  fit  to  the  collected  data  both  are  plotted  on  the  same 
graph,  shown  in  Figure  2. 


Figure  2.  Plot  of  ds/dso  vs.  So,  data  and  fitted  curve  for  a  hypothetical  5  cm  thin  lens. 


10 


3. 


Offset  Measurements 


Actual  complex  lens  systems  have  finite  thickness.  Hence  the  question  arises 
pertaining  to  the  loeations  of  appropriate  origins  from  whieh  to  measure  object  and  image 
distances.  Since  measurements  will  all  be  made  from  the  same  points  with  respeet  to  the 
lens  system  and  image  and  object  planes  there  will  be  no  effect  on  the  differential  values 
Am,  Asi,  and  Aso  used  to  estimate  the  derivatives.  This  constant  offset  between  the  point 
of  measurement  of  objeet  distance  and  the  actual  location  of  the  object  (A)  is  accounted 
for  by  incorporating  it  into  the  function  that  is  used  to  fit  the  data.  The  new  function  is  as 
follows: 

ds,  -f 

dSo 

Where  dsi/dso  is  the  unchanged  y-axis  data,  So  is  the  offset  x-axis  data,  and  A  and/ are  the 
variables  being  ealculated.  A  simulation  of  A  was  incorporated  into  the  thought 
experiment  data  by  artificially  adding  2.4  cm  to  each  simulated  So  data  point.  Results  for 
/  A,  and  the  ealculated  thin  lens  matrix  (Mthm)  closely  mateh  the  expected  values  and  are 
shown  below: 

f  =  5.0012,  A  =  2.3915,  and 


M 


thin 


1.0000  0.0000 

-0.1997  1.0000 


D,  VALIDATION  USING  A  REAL  THIN  LENS 


The  logieal  next  step  in  this  set  of  experiments  is  to  collect  and  analyze  data  on  a 
real  thin  lens  to  ensure  that  the  technique  works  in  practice.  Data  were  collected  and 
analyzed  for  a  thin  glass  lens  whose  nominal  focal  length  was  known  a  priori  to  be  about 
10  em. 

1,  Experimental  Setup 

This  experiment  was  set  up  on  a  standard  optieal  bench.  A  light  source  with  a 
cross  type  pattern  was  plaeed  on  the  right  hand  side  of  the  bench.  The  surface  of  this 
light  source  was  defined  as  the  object  plane.  An  image  was  formed  on  a  small  sheet  of 
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frosted  glass  placed  on  the  left  hand  side  of  the  bench  and  the  lens  was  placed  between 
the  object  and  image  planes.  Measurements  were  made  from  the  center  of  the  lens  to  the 
surfaces  of  the  frosted  glass  and  the  light  source.  Defining  image  and  object  distance  in 
this  way  most  closely  approximates  the  parameters  for  an  ideal  thin  lens. 

2.  Data  and  Results 

The  data  collected  from  this  experiment  is  located  in  appendix  A.  After  the  initial 
calculations  were  performed  it  became  clear  that  some  anomalous  data  points  were 
causing  inaccurate  results.  Spurious  data  points  were  identified  and  eliminated.  All 
derivative  data  whose  values  were  less  than  or  equal  to  -5  were  discarded.  Additionally, 
the  asymptotic  So  data  above  40  cm  did  not  contribute  any  significant  information  and  so 
were  discarded.  Figure  3  is  a  plot  of  all  collected  data  and  graphically  shows  which  data 
points  were  discarded. 


dsi/dso  vs.  Average  Object  Distance  for  a  10  cm  Thin  Lens 


Initially  this  data  set  looked  acceptable  and  returned  credible  results.  The 
resulting  focal  length  was  10.085  cm,  close  to  the  expected  value.  However,  when  the 
average  object  distance  data  were  artificially  offset,  as  in  the  thought  experiment,  the 
calculated  focal  length  was  slightly  greater  than  7  cm,  substantially  different  from  the 
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expected  value.  This  difference  is  due  to  inaccuracies  in  measurement  while  taking  data. 
While  the  actual  distances  were  measured  to  within  a  millimeter,  the  error  associated  with 
light  source  and  lens  positioning  required  for  image  formation  was  much  greater  as  object 
distance  approached  the  lens’  focal  length.  After  the  indicated  data  were  discarded  the 
following  results  were  obtained: 


M,un  = 


1.0000  0.0000 

-0.1021  1.0000 


,  and 


Calculated  Offset  =  -0.9197 


^offset 


1.0000  0.0000 
-0.1021  1.0000 
Calculated  Offset  =  4.0874  cm 


Initial  calculations  on  the  retained  data  were  performed  using  the  function  that 
assumed  no  offset.  Focal  length  yielded  from  the  first  set  of  calculations  was  9.2332  cm. 
While  somewhat  lower  than  expected  this  value  was  accepted  as  a  reasonable  validation 
of  the  technique.  Results  from  the  data  after  a  5  cm  offset  had  been  applied  were  as  seen 
above.  The  calculated  offset  was  about  1  cm  less  than  expected,  and  the  focal  length  was 
approximately  6  mm  greater  than  the  previous  result.  The  most  likely  cause  of  this 
discrepancy  was  related  to  the  assumption  that  the  principal  planes  for  the  thin  lens  were 
located  at  the  center  of  the  lens.  An  ideal  thin  lens  is  assumed  to  be  infinitesimally  thin, 
with  all  of  its  cardinal  points  co-located  on  the  lens  plane.  This  is  obviously  not  the  case 
with  a  real  thin  lens  of  finite  thickness  but  becomes  a  good  assumption  as  image  and 
object  distances  become  large  compared  to  the  focal  length  of  the  lens.  The  calculated 
focal  lengths  agree  to  within  a  millimeter  for  the  two  cases  shown  above,  when  the  offset 
function  is  used  to  fit  both  sets  of  data.  Additionally,  when  the  calculated  offsets  are 
combined  the  artificially  entered  offset  of  5  cm  is  recovered,  as  expected. 

The  thin  lens  thought  experiment  simulations  and  laboratory  data  analysis  were 
deemed  successful  and  a  validation  of  the  technique  for  determining  the  lens  matrix. 
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III.  DETERMINATION  OF  MATRIX  ELEMENTS  FOR 
COMPLEX  OPTICAL  SYSTEMS 


A,  COMPLEX  LENS  SYSTEM  SIMULATIONS 

Similar  to  the  previous  thought  experiments,  a  notional  complex  lens  system  was 
‘designed’  to  verify  that  the  lens  system  matrix  elements  could  be  calculated  from 
measurable  data.  The  first  complex  lens  system  was  comprised  of  three  thin  lenses  with 
constant  separations.  Its  nominal  input  and  output  planes  were  located  at  the  outside 
faces  of  the  individual  lenses  that  comprised  the  system  (see  Figure  4).  A  second 
complex  lens  system  was  considered  which  incorporated  idealized  offsets  at  each  end  of 
the  system.  This  second  hypothetical  system  is  shown  in  Figure  5.  These  mental 
constructs  were  used  to  verily  that  the  matrix  element  determination  method  could  be 
trusted  to  analyze  actual  data  taken  on  the  LINUS  lenses. 


Object  Distance 


t  A 

fl  f2 


Image  Distance 


Y 


Optical  Axis 


Figure  4.  First  notional  complex  lens  system. 
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Figure  5.  Second  notional  complex  lens  system. 

Data  Generation  for  the  Simulations 


Each  of  the  complex  systems  above  has  a  characteristic  matrix.  Using  the 
nomenclature  Mi,  M2,  and  M3  for  the  component  thin  lens  matrices,  Ti,  and  T2  for  the 
lens  separation  matrices,  and  Tqi,  and  T02  for  the  outer  offset  matrices,  we  can  predict 
the  system  matrices  for  each  complex  lens  system  as  follows: 
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fl  ^2  _j_  ^1^2  ^\fl  ^ifl 

fl  fJl 


d-\f2  ~*~  ^2/2  ^1^2 
/2 

^1^2  ~  ^1/2  ~  ^1/2  ~  ^2/2 

Uf, 


a  p\ 

,  ana 

r 


The  expanded  form  of  Msystem2  is  not  shown  beeause  of  its  size.  Hypothetical 
values  defined  for  focal  length,  separation  distance,  and  offset  distance  are  shown  in 
Table  1. 


Table  1.  Hypothetical  values  for  the  notional  systems  shown  in  Figures  4  and  5. 
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The  following  system  matriees  were  calculated  based  on  the  hypothetical  values 
given  in  Table  1. 


M 


system^ 


M 


system2 


0.9322 

1.0063 

-0.0641 

0.9722 

0.8853 

2.7041 

-0.0641 

0.9338 

,  and 


When  these  lens  system  matrices  are  appropriately  combined  with  the  transfer 
matrices  associated  with  the  object  and  image  distances  the  optical  system  matrix  is 
obtained.  The  ‘B’  element  of  this  matrix  is,  by  definition,  equal  to  zero  since  a  finite- 
conjugate  image  is  being  formed.  A  direct  relationship  between  image  and  object 
distance  is  revealed  as  shown  below. 


M 


optical  system 


fl 

■M, 

n 

vO 

lens  system 

vO 

ij 

and 


1.0063  +  0.9722-5, 

0.9622-0.0641-5, 

for  the  lens  system  shown  in  Figure  4. 

A  similar  relationship  follows  for  the  second  lens  system.  Data  were  generated  by 
incrementing  the  simulated  image  distance  from  15.01  to  50  cm  in  0.01  cm  increments 
and  calculating  the  subsequent  object  distances  per  the  formula  above.  Comparable 
image  and  object  distances  were  expected  to  be  measurable  in  the  laboratory  for  various 
wavelengths  of  light.  From  these  measurements  the  change  in  image  distance  with 
respect  to  the  change  in  object  distance  can  be  determined  and  plotted  against  the  average 
object  distance,  as  was  the  case  for  the  thin  lens,  as  described  in  Chapter  II. 

2,  Matrix  Element  Determination 

The  data  generated  in  the  simulations  were  plotted  and  fitted  to  the  same  function 
used  before  for  the  thin  lens  with  a  constant  offset.  Ideally,  this  method  should  return  the 

focal  length  of  the  complex  lens  as  measured  from  the  adjacent  principle  plane  and  the 
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distance  between  the  input  plane  and  that  prineiple  plane  (fi  and  r,  respeetively,  as  shown 
in  Figure  6).  Figure  6  [adapted  from  ref  3,  p.  75]  and  Table  2  [adapted  from  ref  3,  p.  77] 
show  the  geometric  relationships  between  lens  system  matrix  elements,  system  foeal 
lengths  (F),  principle  planes  (H),  and  input  and  output  planes.  The  nodal  points  Ni  and 
N2  are  not  shown. 


Input  Output 

Plane  i  ,  Plane 


Principal  Planes 


Figure  6.  Loeations  of  eardinal  points  Fi,  F2,  Hi,  and  H2  in  a  eomplex  optical  system. 
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Measured  from 
input  and 
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Measured  from 
principal  planes 


Table  2.  Cardinal  point  loeations  defined  in  terms  of  lens  system  matrix  elements. 

Note  that  P  cannot  be  determined  using  these  relationships.  However,  reealling 
that  the  optical  system  and  lens  system  matrix  determinants  must  equal  the  ratio  of  the 
indiees  of  refraction  of  the  media  at  the  system’s  input  and  output,  whieh  is  one  for  this 
experiment,  P  can  be  determined  onee  the  other  matrix  elements  are  known. 
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Results  from  the  data  generated  for  these  notional  lens  systems  were  as  expeeted. 
The  least  squares  curve-fitting  function  returned  -f  and  -r  from  which,  y  and  5  were 
determined.  Since  the  relationship  between  r  and  s  involves  three  unknowns  and  two 
equations  there  is  no  way  to  determine  a  with  only  the  data  at  hand.  It  was  determined 
that  a  second  set  of  data  in  the  reverse  direction  through  the  lens  would  return  the  values 
necessary  for  calculating  a. 


The  ‘reverse’  data  set  was  generated  using  a  relationship  between  image  and 
object  distance  similar  to  the  one  used  while  generating  the  first  set.  The  relationship  was 
obtained  by  determining  a  new  lens  system  matrix  and  then  following  the  same  procedure 
as  previously  discussed.  The  forward  and  reverse  lens  matrices  are  very  closely  related, 
as  shown  below  using  those  from  the  second  notional  system. 


= 


0.8853  2.7041' 

^-0.0641  0.9338y 
^  0.9338  2.7041^ 


and 


^-0.0641  0.8853y 


This  second  data  set  returned  -fi  and  -s,  from  which  a  and  a  second  value  for  y 
were  calculated.  The  second  determination  of  y  provides  an  indication  of  an  appropriate 
level  of  confidence  in  the  measured  data.  If  the  two  values  are  very  close  then  a  high 
level  of  confidence  is  warranted.  After  analysis  was  completed  on  the  generated  data  the 
following  results  were  seen; 


M 


known 


"  0.8853  2.704r 

,-0.0641  0.9338^ 


M 


calculated 


'  0.8852  2.7042^ 

,-0.0641  0.9338^ 


The  consistency  is  excellent. 

B,  LINUS  CAMERA  OBJECTIVE  CHARACTERIZATION  IN  VISIBLE 
LIGHT 

LINUS’  lens  systems  are  aplanats  comprised  of  two  and  three  individual  lenses  of 
unknown  composition  with  unknown  surface  curvatures  and  element  separations  placed 

in  black  barrel  mounts.  Each  barrel  is  marked  only  with  regard  to  the  proper  light 
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propagation  direction.  The  lenses  are  coated  with  an  anti-reflective  material  designed  for 
optimum  throughput  at  near-ultraviolet  wavelengths.  The  camera  objective  for  LINUS 
was  initially  characterized  experimentally  using  four  different  colors  of  visible  light.  In 
these  experiments,  object  distance  was  treated  as  the  independent  variable  and  shortened 
in  0.5  cm  increments  from  a  maximum  value  of  26  cm  down  to  12  cm.  As  can  be  seen  in 
Figure  3,  the  most  critical  data  points  are  those  near  the  ‘knee’  of  the  dsi/dso  vs.  So  curve 
since  they  essentially  define  the  curve.  Data  was  collected  using  red,  yellow,  green,  and 
blue-green  band  pass  filters  and  analyzed  using  the  techniques  previously  discussed. 

During  initial  data  collection  and  analysis  it  was  discovered  that  the  anti-reflective 
lens  coating  did  not  work  at  the  visible  wavelengths.  Lens  flare,  due  to  multiple 
reflections  inside  the  lens  system,  resulted  in  four  images  being  formed  by  the  system 
instead  of  one.  One  image  was  virtual  and  the  other  three  were  real.  The  first  real  image 
was  very  close  to  the  output  plane  and  could  be  disregarded.  The  two  remaining  images 
were  separated  by  approximately  1  cm  and  it  could  not  be  readily  determined  which  was 
the  desired  image  and  which  was  an  artifact  of  the  lens  flare.  Consequently  data  were 
taken  and  analyzed  for  both  sets  of  images.  Subsequent  data  analysis  revealed  that  the 
images  farthest  from  the  output  plane  returned  focal  lengths  as  measured  from  the 
principal  planes  that  were  consistent  while  the  images  closer  to  the  output  plane  returned 
inconsistent  values.  This  fact  indicated  that  the  image  closer  to  the  output  plane  was  an 
artifact  of  lens  flare  and  should  be  neglected. 

1.  Red  Filter  Results 

Data  taken  using  the  red  filter  resulted  in  a  focal  length  from  the  principal  planes 
of  10.98  cm,  a  forward  focal  length  of  9.15  cm,  and  a  reverse  focal  length  of  7.7  cm.  The 
lens  system  matrix  was: 


f  0.8328  4.5682^ 

,-0.0910  0.7014^ 
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2. 


Yellow  Filter  Results 


Data  taken  using  the  yellow  filter  resulted  in  a  foeal  length  from  the  prineipal 
planes  of  10.34  em,  a  forward  foeal  length  of  9.5  em,  and  a  reverse  foeal  length  of  8.17 
em.  The  lens  system  matrix  was: 


M 


yellow 


'  0.9188  2.830r 

^-0.0967  0.7904^ 


3,  Green  Filter  Results 

Data  taken  using  the  green  filter  resulted  in  a  foeal  length  from  the  prineipal 
planes  of  10.37  em,  a  forward  foeal  length  of  9.56  em,  and  a  reverse  foeal  length  of  8.08 
em.  The  lens  system  matrix  was: 


^  0.9226  2.9095^ 

^-0.0965  0.7797^ 


4.  Blue-Green  Filter  Results 

Data  taken  using  the  blue-green  filter  resulted  in  a  foeal  length  from  the  prineipal 
planes  of  10.8  em,  a  forward  foeal  length  of  9.33  em,  and  a  reverse  foeal  length  of  7.77 
em.  The  lens  system  matrix  was: 


f  0.8636  4.09663 

^-0.0926  0.7188^ 
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IV.  DETERMINATION  OF  MATRIX  ELEMENTS  FOR  THE 
LINUS  PRIMARY  AND  CAMERA  OBJECTIVE  LENSES  AT 
ULTRAVIOLET  WAVELENGTHS 


A,  EXPERIMENTAL  SETUP 

This  portion  of  the  thesis  experiment  was  performed  using  a  platinum  hollow 
cathode  lamp  to  illuminate  a  simple  sewing  needle.  The  shadow  image  of  the  needle  was 
then  focused  onto  the  image  intensifier  of  LINUS ’s  UV  camera  through  a  band-pass 
filter.  A  photograph  of  the  experimental  layout  is  shown  in  Figure  7. 


Complex  Lens 
System 


Needle 


Platinum  Hollow 
Cathode  Lamp 


Bending  Mirror 


Lamp  Power 
Supply 


LINUS  Camera 


UV  Filter 


Figure  7.  Physical  layout  used  for  measurements  taken  with  the  camera  objective  lens. 


WinView32,  a  Windows-based  image  acquisition  and  camera  control  program 
provided  by  the  camera  manufacturer,  provided  the  user-camera  interface  to  allow 
focusing  of  the  needle.  This  program  allows  the  user  to  control  the  exposure  time  for  the 
camera,  determine  whether  the  image  is  inverted  or  normal,  and  control  most  of  the  other 
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advanced  features  of  the  eamera.  For  the  purposes  of  this  experiment  the  camera  was  run 
primarily  in  the  ‘focus’  mode  (a  free-running  multiple  frame  data  acquisition  mode)  and 
exposure  time  was  adjusted  between  0.25  and  4  seconds  as  necessary  to  position  the 
needle  correctly  for  image  formation.  A  screen  eapture  of  this  interfaee  with  an  image  of 
the  needle  was  taken  and  is  shown  in  the  figure  below. 


Figure  8.  The  eamera  interface  program  eontrol  window  and  a  needle  image. 

B.  DATA  GENERATION  AND  ANALYSIS 

Image  distanees  were  initially  incremented  in  1  em  steps  to  generate  the 
eorresponding  object  distances.  As  the  change  in  objeet  distance  became  less  than  unity 
the  image  distanee  inerements  were  inereased  to  2,  3,  and  4  em.  This  teehnique  was 
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employed  to  reduee  error  propagation  as  the  data  were  manipulated  into  the  form 
required  for  computer  analysis.  Data  was  taken  in  the  forward  and  reverse  directions  for 
both  the  camera  and  the  primary  objective  lenses  using  four  different  ultraviolet  filters 
centered  at  220,  300,  334,  and  370  nm  with  approximately  10  nm  bandwidths. 

Each  data  set  was  input  to  a  spreadsheet  for  further  manipulation.  Offsets  from 
the  measurement  points  were  accounted  for  and  final  values  for  object  distance  were 
generated.  Next,  values  for  Asi,  Aso,  average  So,  and  Asi/Aso  were  calculated.  All  of  these 
data  are  included  in  Appendix  A.  The  Matlab  code  generated  to  complete  the  analysis  is 
included  in  Appendix  B. 

Initial  analysis  was  performed  using  the  method  described  in  Chapter  III  and  was 
found  to  be  unsatisfactory.  There  was  an  excessive  amount  of  data  scatter  in  the  As/Aso 
vs.  Avg.  So  plots,  giving  somewhat  dubious  results.  In  an  effort  to  reduce  the  amount  of 
uncertainty  associated  with  analyzing  the  data,  So  was  plotted  against  Si  and  fit  to  a 
general  hyperbolic  function  that  would  be  expected  from  theoretical  considerations.  This 
function  was  then  algebraically  manipulated  into  the  correct  form,  differentiated  and  used 
to  generate  values  for  dsi/dso  from  object  distance  data.  The  function  and  its  derivative 
are  shown  below. 


from  which. 


as^  +  b 
cs.  +  d 


ds^  -b 
a -cs^' 


dsi  _  ad  -  be 
dSo  (cs^-af 


Optimal  values  for  a,  b,  c,  and  d  were  calculated  by  the  Matlab  curve  fitting  function 
Tsqcurvefif  and  then  used  with  the  object  distance  data  already  present  to  generate  a 
refined  data  set.  This  data  set  was  then  fit  using  the  function  defined  earlier,  which 
returned  the  focal  lengths  as  measured  from  the  principal  planes  and  the  offset  from  the 
input  and  output  planes.  Values  for  forward  and  reverse  focal  length  as  measured  from 
the  respective  barrel  face  and  the  lens  system  matrix  were  also  calculated  and  returned. 
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C.  RESULTS  FOR  THE  PRIMARY  OBJECTIVE  LENS  SYSTEM 
1,  220  nm  Wavelength 


Data  taken  using  the  220  nm  bandpass  filter  resulted  in  foeal  lengths  from  the 
prineipal  planes  of  18.2  and  24.3  em,  a  forward  foeal  length  of  18.9  em,  and  a  reverse 
foeal  length  of  3 1 .5  em.  The  system  lens  matrix  was; 

^  0.8883  -6.7445^ 

^-0.0471  1.4833  / 


M 

Primary-220 


Plots  of  image  distance  vs.  object  distance  and  ds/dso  vs.  object  distance 
including  the  fitted  curve  are  shown  in  Figure  9. 


Forward  Data 


Reverse  Data 


Forward  Data  (Asymptote  is  reverse  focal  length) 


Reverse  Data  (Asymptote  is  forward  focal  length) 


Figure  9.  Plots  of  Si  vs.  So  and  dsi/dso  vs.  So  for  the  primary  objective  at  220  nm. 
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2. 


300  nm  Wavelength 


Data  taken  using  the  300  nm  bandpass  filter  resulted  in  focal  lengths  from  the 
principal  planes  of  26.2  and  26.5  cm,  a  forward  focal  length  of  23.4  cm,  and  a  reverse 
focal  length  of  25.7  cm.  The  system  lens  matrix  was: 


M 


Primary-300 


0.8903  3.4341 

-0.0380  0.9768 


Plots  of  image  distance  vs.  object  distance  and  ds/dso  vs.  object  distance 
including  the  fitted  curve  are  shown  in  Figure  10. 


Forward  Data 


Reverse  Data 


Forward  Data  (Asymptote  is  reverse  focal  length) 


Reverse  Data  (Asymptote  is  fonward  focal  length) 


Figure  10.  Plots  of  Si  vs.  So  and  dsi/dso  vs.  So  for  the  primary  objective  at  300  nm. 
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3. 


334  nm  Wavelength 


Data  taken  using  the  334  nm  bandpass  filter  resulted  in  foeal  lengths  from  the 
prineipal  planes  of  27.2  and  26.8  cm,  a  forward  focal  length  of  24.1  cm,  and  a  reverse 
focal  length  of  25. 1  cm.  The  system  lens  matrix  was: 


M 


Primary-334 


0.8916  4.6043 

-0.0370  0.9303 


Plots  of  image  distance  vs.  object  distance  and  ds/dso  vs.  object  distance 
including  the  fitted  curve  are  shown  in  Figure  11. 
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4. 


370  nm  Wavelength 


Data  taken  using  the  370  nm  bandpass  filter  resulted  in  foeal  lengths  from  the 
prineipal  planes  of  19.9  and  26.5  cm,  a  forward  focal  length  of  21.4  cm,  and  a  reverse 
focal  length  of  35.9  cm.  The  system  lens  matrix  was: 


M 


Primary-370 


0.9242  -10.0199 

-0.0431  1.5498 


Plots  of  image  distance  vs.  object  distance  and  ds/dso  vs.  object  distance 
including  the  fitted  curve  are  shown  in  Figure  12. 


Forward  Data 


Reverse  Data 


Forward  Data  (Asymptote  is  reverse  focal  length) 


Reverse  Data  (Asymptote  is  forward  focal  length) 


Figure  12.  Plots  of  Si  vs.  Sq  and  dsi/dso  vs.  Sq  for  the  primary  objective  at  370  nm. 
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D.  RESULTS  FOR  THE  CAMERA  OBJECTIVE  LENS  SYSTEM 
1,  220  nm  Wavelength 

Data  taken  using  the  220  nm  bandpass  filter  resulted  in  foeal  lengths  from  the 
prineipal  planes  of  48.5  and  39.8  em,  a  forward  foeal  length  of  58.8  em,  and  a  reverse 
foeal  length  of  41 .8  em.  The  system  lens  matrix  was; 


M 


Camera-220 


M.3319  -11.4763^ 

,-0.0226  0.9460  , 


Plots  of  image  distance  vs.  object  distance  and  ds/dso  vs.  object  distance 
including  the  fitted  curve  are  shown  in  Figure  13. 


Forward  Data 


Forward  Data  (Asymptote  is  reverse  focal  length) 


Reverse  Data  Reverse  Data  (Asymptote  is  forward  focal  length) 


Figure  13.  Plots  of  Si  vs.  Sq  and  ds/dso  vs.  So  for  the  camera  objective  at  220  nm. 


30 


2,  300  nm  Wavelength 

Data  taken  using  the  300  nm  bandpass  filter  resulted  in  foeal  lengths  from  the 
prineipal  planes  of  50.7  and  50.2  cm,  a  forward  focal  length  of  53.0  cm,  and  a  reverse 
focal  length  of  5 1 .6  cm.  The  system  lens  matrix  was: 


M 


Camera-300 


"  1.0493  -3.7035^ 

,-0.0198  1.0230  ^ 


Plots  of  image  distance  vs.  object  distance  and  ds/dso  vs.  object  distance 
including  the  fitted  curve  are  shown  in  Figure  14. 


Forward  Data 


Forward  Data  (Asymptote  is  reverse  focal  length) 


Object  Distance  (cm) 


Reverse  Data 


Reverse  Data  (Asymptote  is  forward  focal  length) 


Figure  14.  Plots  of  Si  vs.  Sq  and  dsi/dso  vs.  So  for  the  camera  objective  at  300  nm. 
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3,  334  nm  Wavelength 

Data  taken  using  the  334  nm  bandpass  filter  resulted  in  foeal  lengths  from  the 
prineipal  planes  of  51.7  and  52.6  cm,  a  forward  focal  length  of  52.0  cm,  and  a  reverse 
focal  length  of  52.9  cm.  The  system  lens  matrix  was: 


M. 


Camera-334 


0.9968  -0.5639 

-0.0192  1.0141 


Plots  of  image  distance  vs.  object  distance  and  ds/dso  vs.  object  distance 
including  the  fitted  curve  are  shown  in  Figure  15. 


Forward  Data 


Forward  Data  (Asymptote  is  reverse  focal  length) 


Object  Distance  (cm) 
Reverse  Data 


Reverse  Data  (Asymptote  is  forward  focal  length) 


Figure  15.  Plots  of  Si  vs.  So  and  dsi/dso  vs.  So  for  the  camera  objective  at  334  nm. 
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4,  370  nm  Wavelength 

Data  taken  using  the  370  nm  bandpass  filter  resulted  in  focal  lengths  from  the 
principal  planes  of  51.9  and  55.2  cm,  a  forward  focal  length  of  55.1  cm,  and  a  reverse 
focal  length  of  57.0  cm.  The  system  lens  matrix  was: 


M 


Camera-370 


^  1.0281  -5.0892^ 

,-0.0187  1.0651  , 


Plots  of  image  distance  vs.  object  distance  and  ds/dso  vs.  object  distance 
including  the  fitted  curve  are  shown  in  Figure  16. 
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Figure  16.  Plots  of  Si  vs.  Sq  and  dsi/dso  vs.  So  for  the  camera  objective  at  370  nm. 
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E,  ERROR  ANALYSIS 

Typically,  error  bars  are  used  when  plotting  data  against  a  curve  or  function  that 
has  been  fitted  to  that  data.  In  this  case  the  error  bars  were  not  plotted  beeause  they  were 
approximately  the  same  size  as  the  eireles  represent  the  data  points  themselves. 
Measurements  eollected  using  the  eamera  objective  are  naturally  larger  due  to  the  longer 
foeal  length  of  that  lens  and  the  fact  that  a  bending  mirror  had  to  be  used  in  order  to 
perform  the  experiment.  As  the  number  of  measurements  inereases  due  to  an  inereased 
number  of  components  in  the  optical  path  the  errors  assoeiated  with  eaeh  measurement 
compound.  The  error  associated  with  eaeh  individual  measurement  of  Si  was  14  of  one 
gradation  of  the  measuring  deviee  used.  Since  that  measurement  was  made  in  one  step 
for  the  primary  objeetive,  that  number  is  the  error  associated  with  those  measurements. 
However,  the  camera  objeetive  required  image  distances  longer  than  the  measuring 
device  being  used.  In  order  to  accomplish  this,  an  optical  mounting  base  was  attached  to 
the  optical  table  at  a  constant  distance  from  the  base  of  the  eamera  and  Si  was  measured 
from  this  point,  effectively  doubling  the  number  of  measurements  required  to  determine 
Si.  The  error  (5)  assoeiated  with  these  values  of  Si  was  determined  using  the  following 
standard  equation: 

^ total  =  +  ^2  ■ 

Finally,  as  the  filter  wavelengths  approaehed  those  of  visible  light  the  image  grew  less 
sharp  when  focused.  This  is  not  unexpected,  as  the  lens  systems  were  designed  to  have 
minimal  aberrations  at  UV  wavelengths,  with  optimum  performance  at  A,~300  nm.  Errors 
assoeiated  with  image  and  objeet  distanee  measurements  are  given  in  the  table  below. 


Lens  System  and  Filter 

5Si 

5so 

Primary  Objective  (220  nm) 

+/-  0.5  mm 

+/-  1  mm 

Primary  Objective  (300  nm) 

+/-  0.5  mm 

+/-  1  mm 

Primary  Objeetive  (334  nm) 

+/-  0.5  mm 

+/-  1  mm 

Primary  Objective  (370  nm) 

+/-  0.5  mm 

+/-  2  mm 
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Camera  Objective  (220  nm) 

+/-  0.71  mm 

+/-  2  mm 

Camera  Objective  (300  nm) 

+/-  0.71  mm 

+/-  2  mm 

Camera  Objective  (334  nm) 

+/-  0.71  mm 

+/-  2  mm 

Camera  Objective  (370  nm) 

+/-  0.71  mm 

+/-  3  mm 

Table  3.  Image  and  object  distance  measurement  errors. 


Recall  that  the  confidence  level  in  the  results  of  the  analysis  is  directly  related  to 
the  focal  lengths  as  measured  from  the  principal  planes.  These  focal  lengths  should,  in 
theory,  be  exactly  equal  in  magnitude  and  as  the  difference  in  their  magnitudes  grows,  in 
reality,  the  level  of  confidence  in  the  results  decreases  proportionally.  The  results  given 
earlier  show  a  high  level  of  confidence  in  the  300  and  334  nm  results  from  both  lens 
systems  and  a  much  lower  level  confidence  in  the  220  and  370  nm  results.  The  fact  that 
the  same  trend  is  seen  from  both  lens  systems  indicates  that  these  less  reliable  results  are 
somehow  linked  to  the  filters  and  does  not  stem  from  the  method,  measurement 
techniques,  lens  systems,  or  other  apparatus  used  to  generate  the  data. 

For  the  purposes  of  this  analysis  the  uncertainty  associated  with  the  forward  and 
reverse  focal  lengths  as  measured  from  the  barrel  is  assumed  to  be  equal  to  the  difference 
in  magnitudes  of  the  calculated  focal  lengths  as  measured  from  the  principal  planes.  The 
forward  and  reverse  focal  lengths  and  their  associated  uncertainties  are  given  in  Tables  4 
and  5. 


Wavelength 

Forward  Focal 

Length 

Reverse  Focal 

Length 

Uncertainty 

220  nm 

18.9  cm 

31.5  cm 

+/-  3.1  cm 

300  nm 

23.4  cm 

25.7  cm 

+/-  0.2  cm 

334  nm 

24.1  cm 

25.1  cm 

+/-  0.5  cm 

370  nm 

21.4  cm 

35.9  cm 

+/-  3.3  cm 

Table  4.  Results  summary  for  the  primary  objective. 
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Wavelength 

Forward 

Focal  Length 

Reverse 

Focal  Length 

Uncertainty 

220  nm 

58.8  cm 

41.8  cm 

+/-  4.9  cm 

300  nm 

53.0  cm 

51.6  cm 

+/-  0.3  cm 

334  nm 

52.0  cm 

52.9  cm 

+/-  0.5  cm 

370  nm 

55.1  cm 

57.0  cm 

+/-  1.7  cm 

Table  5.  Results  summary  for  the  eamera  objeetive. 
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V.  CONCLUSION 


A,  FINAL  RESULTS 

The  primary  goal  of  this  thesis  experiment  was  to  quantify  the  ehromatic 
aberration  charaeteristies  of  the  LINUS  lens  systems.  Due  to  the  lack  of  reliable  results 
from  the  220  and  370  nm  filters  these  results  are  less  accurate  than  originally  hoped.  The 
following  figures  were  generated  by  using  a  simple  spreadsheet;  they  show  the  basic 
trend  of  chromatic  effects  on  focal  length  for  the  primary  and  camera  objective  lenses. 
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B.  RECOMMENDATIONS 


Obviously  the  linear  interpolation  of  two  data  points  is  unambiguous,  but  should 
be  interpreted  with  eaution.  It  is,  however  the  only  option  available  at  this  time.  This 
problem  cannot  be  corrected  until  the  nature  of  the  problems  with  the  220  and  370  nm 
filters  is  understood.  It  is  recommended  that  the  cause  of  errors  in  the  220  and  370  nm 
filter  data  be  explored  and  determined.  Once  the  nature  of  these  errors  is  understood  a 
second  set  of  measurements  should  be  made  and  analyzed.  This  will  allow  for  a  more 
accurate  characterization  of  the  chromatic  aberration  for  these  lenses. 
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APPENDIX  A.  EXPERIMENTAL  DATA 
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Fcm&d 

§ 

^^9 

fimSn 

§ 

in 

§0 

AvqSh 

40 

694 

77.582 

rVa 

rVa 

rVa 

rVa 

40 

58.6 

65.375 

rVa 

rVa 

rVa 

rVa 

41 

56.3 

64.482 

-13.1 

1 

71.032 

-0.07634 

41 

56.1 

62.875 

-2.5 

1 

64.125 

-0.4 

42 

53.4 

61.582 

-2.9 

1 

63.032 

-0.34483 

42 

53.4 

60.175 

-2.7 

1 

61.525 

-0.37037 

43 

49.9 

58.082 

-3.5 

1 

59.832 

-0.28571 

43 

50.9 

57.675 

-2.5 

1 

58.925 

-0.4 

44 

47.9 

56.082 

-2 

1 

57.082 

-0.5 

44 

48.2 

54.975 

-2.7 

1 

56.325 

-0.37037 

45 

46.7 

54.882 

-1.2 

1 

55.482 

-0.83333 

45 

45.9 

52.675 

-2.3 

1 

53.825 

-0.43478 

46 

45.2 

53.382 

-1.5 

1 

54.132 

-0.66667 

46 

44.2 

50.975 

-1.7 

1 

51.825 

-0.58824 

47 

43.9 

52.082 

-1.3 

1 

52.732 

-0.76023 

47 

43.1 

49.875 

-1.1 

1 

50.425 

-0.90909 

48 

41.2 

49.382 

-2.7 

1 

50.732 

-0.37037 

48 

41.7 

48.475 

-1.4 

1 

49.175 

-0.71429 

49 

40.4 

48.582 

-0.8 

1 

48.982 

-1.25 

49 

41.1 

47.875 

-0.6 

1 

48.175 

-1.66667 

50 

39.2 

47.382 

-1.2 

1 

47.982 

-0.83333 

50 

40.2 

46.975 

-0.9 

1 

47.425 

-1.11111 

51 

38.9 

47.082 

-0.3 

1 

47.232 

-3.33333 

51 

39.4 

46.175 

-0.8 

1 

46.575 

-1.25 

52 

38.2 

46.382 

-0.7 

1 

46.732 

-1.42857 

52 

38.5 

45.275 

-0.9 

1 

45.725 

-1.11111 

54 

37 

45.182 

-1.2 

2 

45.782 

-1.66667 

54 

37.1 

43.875 

-1.4 

2 

44.575 

-1.42857 

56 

36.3 

44.482 

-0.7 

2 

44.832 

-2.85714 

56 

36 

42.775 

-1.1 

2 

43.325 

-1.81818 

58 

35.1 

43.282 

-1.2 

2 

43.882 

-1.66667 

58 

34.9 

41.675 

-1.1 

2 

42.225 

-1.81818 

60 

33.7 

41.882 

-1.4 

2 

42.582 

-1.42857 

60 

33.7 

40.475 

-1.2 

2 

41.075 

-1.66667 

62 

32.4 

40.582 

-1.3 

2 

41.232 

-1.53846 

62 

32.4 

39.175 

-1.3 

2 

39.825 

-1.53846 

64 

31.1 

39.282 

-1.3 

2 

39.932 

-1.53846 

64 

31.3 

38.075 

-1.1 

2 

38.625 

-1.81818 

66 

30.5 

38.682 

-0.6 

2 

38.982 

-3.33333 

66 

30 

36.775 

-1.3 

2 

37.425 

-1.53846 

60 

29.2 

37.382 

-1.3 

3 

38.032 

-2.30760 

69 

28.9 

35.675 

-1.1 

3 

36.225 

-2.72727 

72 

28.1 

36.282 

-1.1 

3 

36.832 

-2.72727 

72 

27.6 

34.375 

-1.3 

3 

35.025 

-2.30769 

75 

27.1 

35.282 

-1 

3 

35.782 

-3 

75 

26.9 

33.675 

-0.7 

3 

34.025 

4.28571 

78 

26.5 

34.682 

-0.6 

3 

34.982 

-5 

78 

26.1 

32.875 

-0.8 

3 

33.275 

-3.75 

82 

25.6 

33.782 

-0.9 

4 

34.232 

-4.44444 

82 

25.3 

32.075 

-0.8 

4 

32.475 

-5 

86 

25.3 

33.482 

-0.3 

4 

33.632 

-13.3333 

86 

24.9 

31.675 

-0.4 

4 

31.875 

-10 

90 

24.7 

32.882 

-0.6 

4 

33.182 

-6.66667 

90 

24.7 

31.475 

-0.2 

4 

31.575 

-20 

Data  from  the  Primary  Objective  at  300  nm 


Forward 


Reverse 


n/a 


79.882 


75.332 


72.232 


69.632 


66.782 


64.532 


62.332 


60.382 


58.732 


57.432 


56.032 


54.332 


53.032 


52.382 


51.482 


50.432 


49.682 


48.932 


47.682 


46.432 


45.282 


44.082 


43.282 


42.432 


41.532 


40.782 


40.032 


39.382 


38.632 


37.782 


37.032 


n/a 


-0.22727 


-0.21277 


-0.66667 


-0.27027 


-0.5 


-0.4 


-0.52632 


-0.5 


-0.76923 


-0.76923 


-0.66667 


-0.52632 


-1.42857 


-1.66667 


-0.83333 


-1.11111 


-1.66667 


-1.11111 


-1.25 


-2.22222 


-1.42857 


-2 


-3.33333 


-1.81818 


-2.85714 


-2.5 


-2.85714 


-5 


-3.33333 


-5 


-5.71429 


Si 

s^ 

So 

40 

75.8 

82.575 

41 

71 

77.775 

42 

67.9 

74.675 

43 

64 

70.775 

44 

60.8 

67.575 

45 

58.4 

65.175 

46 

56.1 

62.875 

47 

53.6 

60.375 

48 

51.7 

58.475 

49 

50.1 

56.875 

50 

48.4 

55.175 

51 

47.3 

54.075 

52 

46.2 

52.975 

53 

44.8 

51.575 

54 

43.7 

50.475 

55 

42.8 

49.575 

56 

41.9 

48.675 

57 

41.1 

47.875 

58 

40.7 

47.475 

60 

38.8 

45.575 

62 

37.5 

44.275 

64 

36.1 

42.875 

66 

35.1 

41.875 

68 

34.4 

41.175 

70 

33.6 

40.375 

72 

32.8 

39.575 

74 

32.2 

38.975 

76 

31.4 

38.175 

79 

30.6 

37.375 

82 

29.9 

36.675 

86 

29 

35.775 

90 

28.2 

34.975 

n/a 

n/a 

80.175 

-0.20833 

76.225 

-0.32258 

72.725 

-0.25641 

69.175 

-0.3125 

66.375 

-0.41667 

64.025 

-0.43478 

61.625 

-0.4 

59.425 

-0.52632 

57.675 

-0.625 

56.025 

-0.58824 

54.625 

-0.90909 

53.525 

-0.90909 

52.275 

-0.71429 

51.025 

-0.90909 

50.025 

-1.11111 

49.125 

-1.11111 

48.275 

-1.25 

47.675 

-2.5 

46.525 

-1.05263 

44.925 

-1.53846 

43.575 

-1.42857 

42.375 

-2 

41.525 

-2.85714 

40.775 

-2.5 

39.975 

-2.5 

39.275 

-3.33333 

38.575 

-2.5 

37.775 

-3.75 

37.025 

-4.28571 

36.225 

-4.44444 

35.375 

-5 

4^ 

K) 


Data  from  the  Primary  Objective  at  334  nm 


Forward 

Reverse 

Si 

s!o 

So 

ASn 

ASi 

Si 

s^ 

So 

As^ 

ASi 

40 

78.4 

86.582 

n/a 

n/a 

n/a 

n/a 

40 

82 

88.775 

n/a 

n/a 

n/a 

n/a 

41 

74.8 

82.982 

-3.6 

1 

84.782 

-0.27778 

41 

76 

82.775 

-6 

1 

85.775 

-0.16667 

42 

71.1 

79.282 

-3.7 

1 

81.132 

-0.27027 

42 

71.4 

78.175 

-4.6 

1 

80.475 

-0.21739 

43 

67.2 

75.382 

-3.9 

1 

77.332 

-0.25641 

43 

67.2 

73.975 

-4.2 

1 

76.075 

-0.2381 

44 

64 

72.182 

-3.2 

1 

73.782 

-0.3125 

44 

64 

70.775 

-3.2 

1 

72.375 

-0.3125 

45 

61.2 

69.382 

-2.8 

1 

70.782 

-0.35714 

45 

61.1 

67.875 

-2.9 

1 

69.325 

-0.34483 

46 

59 

67.182 

-2.2 

1 

68.282 

-0.45455 

46 

58.9 

65.675 

-2.2 

1 

66.775 

-0.45455 

47 

56.3 

64.482 

-2.7 

1 

65.832 

-0.37037 

47 

56.5 

63.275 

-2.4 

1 

64.475 

-0.41667 

48 

53.9 

62.082 

-2.4 

1 

63.282 

-0.41667 

48 

54 

60.775 

-2.5 

1 

62.025 

-0.4 

49 

52.5 

60.682 

-1.4 

1 

61.382 

-0.71429 

49 

52.1 

58.875 

-1.9 

1 

59.825 

-0.52632 

50 

50.8 

58.982 

-1.7 

1 

59.832 

-0.58824 

50 

50.8 

57.575 

-1.3 

1 

58.225 

-0.76923 

51 

48.9 

57.082 

-1.9 

1 

58.032 

-0.52632 

51 

49.2 

55.975 

-1.6 

1 

56.775 

-0.625 

52 

48 

56.182 

-0.9 

1 

56.632 

-1.11111 

52 

47.9 

54.675 

-1.3 

1 

55.325 

-0.76923 

53 

46.8 

54.982 

-1.2 

1 

55.582 

-0.83333 

53 

46.6 

53.375 

-1.3 

1 

54.025 

-0.76923 

54 

45.2 

53.382 

-1.6 

1 

54.182 

-0.625 

54 

45.5 

52.275 

-1.1 

1 

52.825 

-0.90909 

55 

44.5 

52.682 

-0.7 

1 

53.032 

-1.42857 

55 

44.2 

50.975 

-1.3 

1 

51.625 

-0.76923 

56 

43.8 

51.982 

-0.7 

1 

52.332 

-1.42857 

56 

43.5 

50.275 

-0.7 

1 

50.625 

-1.42857 

57 

42.3 

50.482 

-1.5 

1 

51.232 

-0.66667 

57 

42.6 

49.375 

-0.9 

1 

49.825 

-1.11111 

58 

41.7 

49.882 

-0.6 

1 

50.182 

-1.66667 

58 

41.6 

48.375 

-1 

1 

48.875 

-1 

60 

40.1 

48.282 

-1.6 

2 

49.082 

-1.25 

60 

40.1 

46.875 

-1.5 

2 

47.625 

-1.33333 

62 

38.8 

46.982 

-1.3 

2 

47.632 

-1.53846 

62 

38.8 

45.575 

-1.3 

2 

46.225 

-1.53846 

64 

37.7 

45.882 

-1.1 

2 

46.432 

-1.81818 

64 

37.7 

44.475 

-1.1 

2 

45.025 

-1.81818 

66 

36.6 

44.782 

-1.1 

2 

45.332 

-1.81818 

66 

36.6 

43.375 

-1.1 

2 

43.925 

-1.81818 

68 

35.8 

43.982 

-0.8 

2 

44.382 

-2.5 

68 

35.8 

42.575 

-0.8 

2 

42.975 

-2.5 

70 

34.9 

43.082 

-0.9 

2 

43.532 

-2.22222 

70 

34.5 

41.275 

-1.3 

2 

41.925 

-1.53846 

72 

33.9 

42.082 

-1 

2 

42.582 

-2 

72 

34 

40.775 

-0.5 

2 

41.025 

-4 

74 

33.3 

41.482 

-0.6 

2 

41.782 

-3.33333 

74 

33.3 

40.075 

-0.7 

2 

40.425 

-2.85714 

76 

32.7 

40.882 

-0.6 

2 

41.182 

-3.33333 

76 

32.6 

39.375 

-0.7 

2 

39.725 

-2.85714 

79 

31.9 

40.082 

-0.8 

3 

40.482 

-3.75 

79 

31.6 

38.375 

-1 

3 

38.875 

-3 

82 

31 

39.182 

-0.9 

3 

39.632 

-3.33333 

82 

30.8 

37.575 

-0.8 

3 

37.975 

-3.75 

86 

29.9 

38.082 

-1.1 

4 

38.632 

-3.63636 

86 

29.9 

36.675 

-0.9 

4 

37.125 

-4.44444 

90 

29.2 

37.382 

-0.7 

4 

37.732 

-5.71429 

90 

29 

35.775 

-0.9 

4 

36.225 

-4.44444 

Forward 


Data  from  the  Primary  Objective  at  370  nm 


Si 

As^ 

ASi 

Si 

s^ 

40 

81.8 

89.982 

n/a 

n/a 

n/a 

n/a 

40 

84.2 

90.975 

41 

75.2 

83.382 

-6.6 

1 

86.682 

-0.15152 

41 

78.8 

85.575 

42 

72.2 

80.382 

-3 

1 

81.882 

-0.33333 

42 

75.5 

82.275 

43 

67 

75.182 

-5.2 

1 

77.782 

-0.19231 

43 

71.7 

78.475 

44 

62.1 

70.282 

-4.9 

1 

72.732 

-0.20408 

44 

66.5 

73.275 

45 

59.3 

67.482 

-2.8 

1 

68.882 

-0.35714 

45 

62.5 

69.275 

46 

55.3 

63.482 

-4 

1 

65.482 

-0.25 

46 

59.7 

66.475 

47 

53.4 

61.582 

-1.9 

1 

62.532 

-0.52632 

47 

57.9 

64.675 

48 

51.5 

59.682 

-1.9 

1 

60.632 

-0.52632 

48 

55.9 

62.675 

49 

49.3 

57.482 

-2.2 

1 

58.582 

-0.45455 

49 

53.5 

60.275 

50 

47.8 

55.982 

-1.5 

1 

56.732 

-0.66667 

50 

51.8 

58.575 

51 

47 

55.182 

-0.8 

1 

55.582 

-1.25 

51 

50.4 

57.175 

52 

46 

54.182 

-1 

1 

54.682 

-1 

52 

48.8 

55.575 

53 

44.9 

53.082 

-1.1 

1 

53.632 

-0.90909 

53 

47.6 

54.375 

54 

44.2 

52.382 

-0.7 

1 

52.732 

-1.42857 

54 

46.3 

53.075 

55 

42.8 

50.982 

-1.4 

1 

51.682 

-0.71429 

55 

45.3 

52.075 

56 

42 

50.182 

-0.8 

1 

50.582 

-1.25 

56 

44.7 

51.475 

57 

41.4 

49.582 

-0.6 

1 

49.882 

-1.66667 

57 

43.9 

50.675 

58 

40.6 

48.782 

-0.8 

1 

49.182 

-1.25 

58 

43 

49.775 

59 

40 

48.182 

-0.6 

1 

48.482 

-1.66667 

60 

41.2 

47.975 

60 

39.4 

47.582 

-0.6 

1 

47.882 

-1.66667 

62 

40.1 

46.875 

61 

38.6 

46.782 

-0.8 

1 

47.182 

-1.25 

64 

38.4 

45.175 

62 

37.9 

46.082 

-0.7 

1 

46.432 

-1.42857 

66 

37.1 

43.875 

63 

37.2 

45.382 

-0.7 

1 

45.732 

-1.42857 

68 

36.4 

43.175 

69 

36 

44.182 

-1.2 

6 

44.782 

-5 

70 

35.4 

42.175 

71 

35.2 

43.382 

-0.8 

2 

43.782 

-2.5 

72 

34.5 

41.275 

73 

34.3 

42.482 

-0.9 

2 

42.932 

-2.22222 

74 

33.7 

40.475 

75 

33.6 

41.782 

-0.7 

2 

42.132 

-2.85714 

76 

33.1 

39.875 

77 

32.8 

40.982 

-0.8 

2 

41.382 

-2.5 

79 

32 

38.775 

80 

32 

40.182 

-0.8 

3 

40.582 

-3.75 

82 

31.3 

38.075 

83 

31.3 

39.482 

-0.7 

3 

39.832 

-4.28571 

86 

30.3 

37.075 

86 

30.5 

38.682 

-0.8 

3 

39.082 

-3.75 

90 

29.4 

36.175 

n/a 

n/a 

88.275 

-0.18519 

83.925 

-0.30303 

80.375 

-0.26316 

75.875 

-0.19231 

71.275 

-0.25 

67.875 

-0.35714 

65.575 

-0.55556 

63.675 

-0.5 

61.475 

-0.41667 

59.425 

-0.58824 

57.875 

-0.71429 

56.375 

-0.625 

54.975 

-0.83333 

53.725 

-0.76923 

52.575 

-1 

51.775 

-1.66667 

51.075 

-1.25 

50.225 

-1.11111 

48.875 

-1.11111 

47.425 

-1.81818 

46.025 

-1.17647 

44.525 

-1.53846 

43.525 

-2.85714 

42.675 

-2 

41.725 

-2.22222 

40.875 

-2.5 

40.175 

-3.33333 

39.325 

-2.72727 

38.425 

-4.28571 

37.575 

36.625 

-4.44444 

Forward 


Data  from  the  Camera  Objective  at  220  nm 


Reverse 


Si 

dnrr, 

S_Q 

Sq 

ASn 

ASi 

Si 

dnm 

Sq 

ASn 

ASi 

80 

44.2 

120.9 

127.784 

n/a 

n/a 

n/a 

n/a 

80 

44.5 

121.2 

129.257 

n/a 

n/a 

n/a 

n/a 

81 

41.1 

116.8 

123.684 

■4  A 

1 

125.734 

-0.2439024 

81 

43 

118.7 

126.757 

-2.5 

1 

128.007 

-0.4 

82 

40.5 

115.2 

122.084 

-1.6 

1 

122.884 

-0.625 

82 

40.4 

115.1 

123.157 

-3.6 

1 

124.957 

-0.2777778 

83 

38.4 

112.1 

118.984 

-3.1 

1 

120.534 

-0.3225806 

83 

39.7 

113.4 

121.457 

-1.7 

1 

122.307 

-0.5882353 

84 

36.9 

109.6 

116.484 

-2.5 

1 

117.734 

-0.4 

84 

38 

110.7 

118.757 

-2.7 

1 

120.107 

-0.3703704 

85 

36.2 

107.9 

114.784 

-1.7 

1 

115.634 

-0.5882353 

85 

36.6 

108.3 

116.357 

-2.4 

1 

117.557 

-0.4166667 

86 

35 

105.7 

112.584 

-2.2 

1 

113.684 

-0.4545455 

86 

36.7 

107.4 

115.457 

-0.9 

1 

115.907 

-1.1111111 

87 

34.4 

104.1 

110.984 

-1.6 

1 

111.784 

-0.625 

87 

34 

103.7 

111.757 

-3.7 

1 

113.607 

-0.2702703 

88 

34.2 

102.9 

109.784 

-1.2 

1 

110.384 

-0.8333333 

88 

32.9 

101.6 

109.657 

-2.1 

1 

110.707 

-0.4761905 

89 

33.8 

101.5 

108.384 

-1.4 

1 

109.084 

-0.7142857 

89 

31.9 

99.6 

107.657 

-2 

1 

108.657 

-0.5 

90 

32.4 

99.1 

105.984 

-2.4 

1 

107.184 

-0.4166667 

90 

30 

96.7 

104.757 

-2.9 

1 

106.207 

-0.3448276 

91 

31.8 

97.5 

104.384 

-1.6 

1 

105.184 

-0.625 

91 

29.1 

94.8 

102.857 

-1.9 

1 

103.807 

-0.5263158 

92 

31.8 

96.5 

103.384 

-1 

1 

103.884 

-1 

92 

28.9 

93.6 

101.657 

-1.2 

1 

102.257 

-0.8333333 

93 

31.8 

95.5 

102.384 

-1 

1 

102.884 

-1 

93 

28.3 

92 

100.057 

-1.6 

1 

100.857 

-0.625 

94 

31.9 

94.6 

101.484 

-0.9 

1 

101.934 

-1.1111111 

94 

28.1 

90.8 

98.857 

-1.2 

1 

99.457 

-0.8333333 

95 

31.6 

93.3 

100.184 

-1.3 

1 

100.834 

-0.7692308 

95 

27.8 

89.5 

97.557 

-1.3 

1 

98.207 

-0.7692308 

96 

31.2 

91.9 

98.784 

-1.4 

1 

99.484 

-0.7142857 

96 

27.6 

88.3 

96.357 

-1.2 

1 

96.957 

-0.8333333 

97 

31 

90.7 

97.584 

-1.2 

1 

98.184 

-0.8333333 

97 

27.3 

87 

95.057 

-1.3 

1 

95.707 

-0.7692308 

98 

30.3 

89 

95.884 

-1.7 

1 

96.734 

-0.5882353 

98 

27.2 

85.9 

93.957 

-1.1 

1 

94.507 

-0.9090909 

99 

30.2 

87.9 

94.784 

-1.1 

1 

95.334 

-0.9090909 

99 

28.4 

86.1 

94.157 

0.2 

1 

94.057 

5 

100 

29.6 

86.3 

93.184 

-1.6 

1 

93.984 

-0.625 

100 

29.1 

85.8 

93.857 

-0.3 

1 

94.007 

-3.3333333 

102 

31.1 

85.8 

92.684 

-0.5 

2 

92.934 

-4 

102 

28.4 

83.1 

91.157 

-2.7 

2 

92.507 

-0.7407407 

104 

31.8 

84.5 

91.384 

-1.3 

2 

92.034 

-1.5384615 

104 

29.3 

82 

90.057 

-1.1 

2 

90.607 

-1.8181818 

106 

32 

82.7 

89.584 

-1.8 

2 

90.484 

-1.1111111 

106 

30 

80.7 

88.757 

-1.3 

2 

89.407 

-1.5384615 

108 

31.8 

80.5 

87.384 

-2.2 

2 

88.484 

-0.9090909 

108 

31.9 

80.6 

88.657 

-0.1 

2 

88.707 

-20 

110 

31.9 

78.6 

85.484 

-1.9 

2 

86.434 

-1.0526316 

110 

31.1 

77.8 

85.857 

-2.8 

2 

87.257 

-0.7142857 

112 

32.2 

76.9 

83.784 

-1.7 

2 

84.634 

-1.1764706 

112 

32.2 

76.9 

84.957 

-0.9 

2 

85.407 

-2.2222222 

114 

32.3 

75 

81.884 

-1.9 

2 

82.834 

-1.0526316 

114 

32.7 

75.4 

83.457 

-1.5 

2 

84.207 

-1.3333333 

116 

32.8 

73.5 

80.384 

-1.5 

2 

81.134 

-1.3333333 

116 

35.4 

76.1 

84.157 

0.7 

2 

83.807 

2.85714286 

118 

35.6 

74.3 

81.184 

0.8 

2 

80.784 

2.5 

118 

36.2 

74.9 

82.957 

-1.2 

2 

83.557 

-1.6666667 

120 

37 

73.7 

80.584 

-0.6 

2 

80.884 

-3.3333333 

120 

36.5 

73.2 

81.257 

-1.7 

2 

82.107 

-1.1764706 

122 

37.9 

72.6 

79.484 

-1.1 

2 

80.034 

-1.8181818 

122 

37.1 

71.8 

79.857 

-1.4 

2 

80.557 

-1.4285714 

124 

38.5 

71.2 

78.084 

-1.4 

2 

78.784 

-1.4285714 

124 

38.2 

70.9 

78.957 

-0.9 

2 

79.407 

-2.2222222 

126 

39.6 

70.3 

77.184 

-0.9 

2 

77.634 

-2.2222222 

126 

39.3 

70 

78.057 

-0.9 

2 

78.507 

-2.2222222 

128 

41.1 

69.8 

76.684 

-0.5 

2 

76.934 

-4 

128 

40.4 

69.1 

77.157 

-0.9 

2 

77.607 

-2.2222222 

130 

42.6 

69.3 

76.184 

-0.5 

2 

76.434 

-4 

130 

41.9 

68.6 

76.657 

-0.5 

2 

76.907 

-4 

132 

44.1 

68.8 

75.684 

-0.5 

2 

75.934 

-4 

132 

43 

67.7 

75.757 

-0.9 

2 

76.207 

-2.2222222 

134 

44.9 

67.6 

74.484 

-1.2 

2 

75.084 

-1.6666667 

134 

44.4 

67.1 

75.157 

-0.6 

2 

75.457 

-3.3333333 

136 

46.1 

66.8 

73.684 

-0.8 

2 

74.084 

-2.5 

136 

45.9 

66.6 

74.657 

-0.5 

2 

74.907 

-4 

138 

47.6 

66.3 

73.184 

-0.5 

2 

73.434 

-4 

138 

47.3 

66 

74.057 

-0.6 

2 

74.357 

-3.3333333 

140 

48.9 

65.6 

72.484 

-0.7 

2 

72.834 

-2.8571429 

140 

49.1 

65.8 

73.857 

-0.2 

2 

73.957 

-10 

144 

52.2 

64.9 

71.784 

-0.7 

4 

72.134 

-5.7142857 

144 

51.9 

64.6 

72.657 

-1.2 

4 

73.257 

-3.3333333 

148 

55 

63.7 

70.584 

-1.2 

4 

71.184 

-3.3333333 

148 

55.2 

63.9 

71.957 

-0.7 

4 

72.307 

-5.7142857 

Data  from  the 


Forward 

Si 

dnm 

So 

ASn 

ASi 

ASi/ASn 

80 

72.2 

148.9 

155.784 

n/a 

n/a 

n/a 

n/a 

81 

69.4 

145.1 

151.984 

-3.8 

1 

153.884 

-0.2631579 

82 

66 

140.7 

147.584 

-4.4 

1 

149.784 

-0.2272727 

83 

63.9 

137.6 

144.484 

-3.1 

1 

146.034 

-0.3225806 

84 

61.2 

133.9 

140.784 

-3.7 

1 

142.634 

-0.2702703 

85 

60 

131.7 

138.584 

-2.2 

1 

139.684 

-0.4545455 

86 

57.9 

128.6 

135.484 

-3.1 

1 

137.034 

-0.3225806 

87 

55.9 

125.6 

132.484 

-3 

1 

133.984 

-0.3333333 

88 

54.6 

123.3 

130.184 

-2.3 

1 

131.334 

-0.4347826 

89 

53.1 

120.8 

127.684 

-2.5 

1 

128.934 

-0.4 

90 

51.6 

118.3 

125.184 

-2.5 

1 

126.434 

-0.4 

91 

51.4 

117.1 

123.984 

-1.2 

1 

124.584 

-0.8333333 

92 

50.1 

114.8 

121.684 

-2.3 

1 

122.834 

-0.4347826 

93 

49.2 

112.9 

119.784 

-1.9 

1 

120.734 

-0.5263158 

94 

48.3 

111 

117.884 

-1.9 

1 

118.834 

-0.5263158 

95 

47.5 

109.2 

116.084 

-1.8 

1 

116.984 

-0.5555556 

96 

47.1 

107.8 

1 14.684 

-1.4 

1 

115.384 

-0.7142857 

97 

47.2 

106.9 

113.784 

-0.9 

1 

114.234 

-1.1111111 

98 

46.8 

105.5 

112.384 

-1.4 

1 

113.084 

-0.7142857 

99 

46.2 

103.9 

110.784 

-1.6 

1 

111.584 

-0.625 

100 

45.5 

102.2 

109.084 

-1.7 

1 

109.934 

-0.5882353 

102 

45.4 

100.1 

106.984 

-2.1 

2 

108.034 

-0.952381 

104 

45.2 

97.9 

104.784 

-2.2 

2 

105.884 

-0.9090909 

106 

45 

95.7 

102.584 

-2.2 

2 

103.684 

-0.9090909 

108 

45.2 

93.9 

100.784 

-1.8 

2 

101.684 

-1.1111111 

110 

45.9 

92.6 

99.484 

-1.3 

2 

100.134 

-1.5384615 

112 

45.4 

90.1 

96.984 

-2.5 

2 

98.234 

-0.8 

114 

45.9 

88.6 

95.484 

-1.5 

2 

96.234 

-1.3333333 

116 

46.7 

87.4 

94.284 

-1.2 

2 

94.884 

-1.6666667 

118 

47.1 

85.8 

92.684 

-1.6 

2 

93.484 

-1.25 

o 

CM 

47.6 

84.3 

91.184 

-1.5 

2 

91.934 

-1.3333333 

CM 

CM 

49 

83.7 

90.584 

-0.6 

2 

90.884 

-3.3333333 

CM 

50.4 

83.1 

89.984 

-0.6 

2 

90.284 

-3.3333333 

CO 

CM 

51.1 

81.8 

88.684 

-1.3 

2 

89.334 

-1.5384615 

1  128 

52.1 

80.8 

87.684 

-1 

2 

88.184 

-2 

CO 

o 

52.8 

79.5 

86.384 

-1.3 

2 

87.034 

-1.5384615 

132 

53.7 

78.4 

85.284 

-1.1 

2 

85.834 

-1.8181818 

134 

54.9 

77.6 

84.484 

-0.8 

2 

84.884 

-2.5 

136 

56.2 

76.9 

83.784 

-0.7 

2 

84.134 

-2.8571429 

138 

57.4 

76.1 

82.984 

-0.8 

2 

83.384 

-2.5 

140 

58.5 

75.2 

82.084 

-0.9 

2 

82.534 

-2.2222222 

144 

60.8 

73.5 

80.384 

-1.7 

4 

81.234 

-2.3529412 

148 

63.9 

72.6 

79.484 

-0.9 

4 

79.934 

^.4444444 

Objective  at  300  nm 


Reverse 


Si 

dnm 

Sq 

ASn 

ASi 

ASi/ASn 

80 

74.9 

151.6 

159.657 

n/a 

n/a 

n/a 

n/a 

81 

71.8 

147.5 

155.557 

-4.1 

1 

157.607 

-0.2439024 

82 

68.3 

143 

151.057 

-4.5 

1 

153.307 

-0.2222222 

83 

64.8 

138.5 

146.557 

-4.5 

1 

148.807 

-0.2222222 

84 

63.7 

136.4 

144.457 

-2.1 

1 

145.507 

-0.4761905 

85 

61.1 

132.8 

140.857 

-3.6 

1 

142.657 

-0.2777778 

86 

59.8 

130.5 

138.557 

-2.3 

1 

139.707 

-0.4347826 

87 

57.7 

127.4 

135.457 

-3.1 

1 

137.007 

-0.3225806 

88 

56.6 

125.3 

133.357 

-2.1 

1 

134.407 

-0.4761905 

89 

54.7 

122.4 

130.457 

-2.9 

1 

131.907 

-0.3448276 

90 

53 

119.7 

127.757 

-2.7 

1 

129.107 

-0.3703704 

91 

51 

116.7 

124.757 

-3 

1 

126.257 

-0.3333333 

92 

50.2 

114.9 

122.957 

-1.8 

1 

123.857 

BilA-Cm-fil 

93 

49.5 

113.2 

121.257 

-1.7 

1 

122.107 

-0.5882353 

94 

48.2 

110.9 

118.957 

-2.3 

1 

120.107 

-0.4347826 

95 

47.6 

109.3 

117.357 

-1.6 

1 

118.157 

-0.625 

96 

47.2 

107.9 

115.957 

-1.4 

1 

116.657 

-0.7142857 

97 

47.3 

107 

115.057 

-0.9 

1 

115.507 

-1.1111111 

98 

46.8 

105.5 

113.557 

-1.5 

1 

114.307 

-0.6666667 

99 

46.5 

104.2 

112.257 

-1.3 

1 

112.907 

-0.7692308 

100 

46.3 

103 

111.057 

-1.2 

1 

111.657 

-0.8333333 

CM 

O 

45.8 

100.5 

108.557 

-2.5 

2 

109.807 

-0.8 

104 

45.5 

98.2 

106.257 

-2.3 

2 

107.407 

-0.8695652 

106 

44.8 

95.5 

103.557 

-2.7 

2 

104.907 

-0.7407407 

CO 

o 

44.9 

93.6 

101.657 

-1.9 

2 

102.607 

-1.0526316 

110 

44.4 

91.1 

99.157 

-2.5 

2 

100.407 

-0.8 

112 

45.7 

90.4 

98.457 

-0.7 

2 

98.807 

-2.8571429 

114 

45.9 

88.6 

96.657 

-1.8 

2 

97.557 

-1.1111111 

116 

46 

86.7 

94.757 

-1.9 

2 

95.707 

-1.0526316 

118 

46.8 

85.5 

93.557 

-1.2 

2 

94.157 

-1.6666667 

o 

CM 

47.8 

84.5 

92.557 

-1 

2 

93.057 

-2 

CM 

CM 

48.7 

83.4 

91.457 

-1.1 

2 

92.007 

-1.8181818 

CM 

49.4 

82.1 

90.157 

-1.3 

2 

90.807 

-1.5384615 

CD 

CM 

50.1 

80.8 

88.857 

-1.3 

2 

89.507 

-1.5384615 

CO 

CM 

51.4 

80.1 

88.157 

-0.7 

2 

88.507 

-2.8571429 

CO 

o 

52.3 

79 

87.057 

-1.1 

2 

87.607 

-1.8181818 

132 

53 

77.7 

85.757 

-1.3 

2 

86.407 

-1.5384615 

134 

54.4 

77.1 

85.157 

-0.6 

2 

85.457 

-3.3333333 

136 

55.4 

76.1 

84.157 

-1 

2 

84.657 

-2 

CO 

CO 

56.6 

75.3 

83.357 

-0.8 

2 

83.757 

-2.5 

140 

57.8 

74.5 

82.557 

-0.8 

2 

82.957 

-2.5 

144 

61 

73.7 

81.757 

-0.8 

4 

82.157 

-5 

CO 

63.2 

71.9 

79.957 

-1.8 

4 

80.857 

-2.2222222 

Data  from  the  Camera  Objective  at  334  nm 


Forward 

Reverse 

Si 

^nm 

Slii 

So 

ASr, 

ASi 

Si 

^nm 

s:^ 

So 

ASq 

ASi 

80 

80.2 

156.9 

163.784 

n/a 

n/a 

n/a 

1  1 

80 

83 

159.7 

167.757 

n/a 

n/a 

n/a 

n/a 

81 

78.6 

154.3 

161.184 

-2.6 

1 

162.484 

81 

79 

154.7 

162.757 

-5 

1 

165.257 

-0.2 

82 

74.3 

149 

155.884 

-5.3 

1 

158.534 

82 

75.1 

149.8 

157.857 

-4.9 

1 

160.307 

-0.2040816 

83 

71.5 

145.2 

152.084 

-3.8 

1 

153.984 

83 

73.2 

146.9 

154.957 

-2.9 

1 

156.407 

-0.3448276 

84 

69.4 

142.1 

148.984 

-3.1 

1 

150.534 

84 

70.7 

143.4 

151.457 

-3.5 

1 

153.207 

-0.2857143 

85 

66.6 

138.3 

145.184 

-3.8 

1 

147.084 

85 

67.1 

138.8 

146.857 

-4.6 

1 

149.157 

-0.2173913 

86 

65 

135.7 

142.584 

-2.6 

1 

143.884 

86 

64.1 

134.8 

142.857 

-4 

1 

144.857 

-0.25 

87 

62.5 

132.2 

139.084 

-3.5 

1 

140.834 

-0.2857143 

87 

62.4 

132.1 

140.157 

-2.7 

1 

141.507 

-0.3703704 

88 

60.9 

129.6 

136.484 

-2.6 

1 

137.784 

-0.3846154 

88 

61 

129.7 

137.757 

-2.4 

1 

138.957 

-0.4166667 

89 

59.9 

127.6 

134.484 

-2 

1 

135.484 

-0.5 

89 

59.3 

127 

135.057 

-2.7 

1 

136.407 

-0.3703704 

90 

57.9 

124.6 

131.484 

-3 

1 

132.984 

90 

57.4 

124.1 

132.157 

-2.9 

1 

133.607 

-0.3448276 

91 

56.1 

121.8 

128.684 

-2.8 

1 

130.084 

91 

56.2 

121.9 

129.957 

-2.2 

1 

131.057 

-0.4545455 

92 

54.7 

119.4 

126.284 

-2.4 

1 

127.484 

92 

55.4 

120.1 

128.157 

-1.8 

1 

129.057 

-0.5555556 

93 

54 

117.7 

124.584 

-1.7 

1 

125.434 

93 

55.5 

119.2 

127.257 

1 

127.707 

-1.1111111 

94 

53.2 

115.9 

122.784 

-1.8 

1 

123.684 

94 

54.8 

117.5 

125.557 

-1.7 

1 

126.407 

-0.5882353 

95 

52.3 

114 

120.884 

-1.9 

1 

121.834 

95 

53.6 

115.3 

123.357 

-2.2 

1 

124.457 

-0.4545455 

96 

51.6 

112.3 

119.184 

-1.7 

1 

120.034 

96 

53.3 

114 

122.057 

-1.3 

1 

122.707 

-0.7692308 

97 

51.3 

111 

117.884 

-1.3 

1 

118.534 

-0.7692308 

97 

52.9 

112.6 

120.657 

-1.4 

1 

121.357 

-0.7142857 

98 

50.5 

109.2 

116.084 

-1.8 

1 

116.984 

-0.5555556 

98 

52.1 

110.8 

118.857 

-1.8 

1 

119.757 

-0.5555556 

99 

49.8 

107.5 

114.384 

-1.7 

1 

115.234 

-0.5882353 

99 

51.9 

109.6 

117.657 

-1.2 

1 

118.257 

-0.8333333 

100 

49.5 

106.2 

113.084 

-1.3 

1 

113.734 

100 

51.7 

108.4 

116.457 

-1.2 

1 

117.057 

-0.8333333 

102 

48.8 

103.5 

110.384 

-2.7 

2 

111.734 

WikBirmA 

102 

49.9 

104.6 

112.657 

-3.8 

2 

114.557 

-0.5263158 

104 

48.6 

101.3 

108.184 

-2.2 

2 

109.284 

104 

48.9 

101.6 

109.657 

-3 

2 

111.157 

-0.6666667 

106 

48.2 

98.9 

105.784 

-2.4 

2 

106.984 

Bmi^ 

106 

48.2 

98.9 

106.957 

-2.7 

2 

108.307 

-0.7407407 

108 

49 

97.7 

104.584 

-1.2 

2 

105.184 

-1.6666667 

108 

47.8 

96.5 

104.557 

-2.4 

2 

105.757 

-0.8333333 

110 

48.7 

95.4 

102.284 

-2.3 

2 

103.434 

-0.8695652 

110 

48.7 

95.4 

103.457 

-1.1 

2 

104.007 

-1.8181818 

112 

49 

93.7 

100.584 

-1.7 

2 

101.434 

-1.1764706 

112 

49 

93.7 

101.757 

-1.7 

2 

102.607 

-1.1764706 

114 

49.2 

91.9 

98.784 

-1.8 

2 

99.684 

-1.1111111 

114 

49.4 

92.1 

100.157 

-1.6 

2 

100.957 

-1.25 

116 

49.9 

90.6 

97.484 

-1.3 

2 

98.134 

-1.5384615 

116 

50 

90.7 

98.757 

-1.4 

2 

99.457 

-1.4285714 

118 

50.6 

89.3 

96.184 

-1.3 

2 

96.834 

-1.5384615 

118 

50.7 

89.4 

97.457 

-1.3 

2 

98.107 

-1.5384615 

120 

51.4 

88.1 

94.984 

-1.2 

2 

95.584 

-1.6666667 

120 

51.2 

87.9 

95.957 

-1.5 

2 

96.707 

-1.3333333 

122 

51.9 

86.6 

93.484 

-1.5 

2 

94.234 

-1.3333333 

122 

51.9 

86.6 

94.657 

-1.3 

2 

95.307 

-1.5384615 

124 

52.6 

85.3 

92.184 

-1.3 

2 

92.834 

-1.5384615 

124 

52.6 

85.3 

93.357 

-1.3 

2 

94.007 

-1.5384615 

126 

53.3 

84 

90.884 

-1.3 

2 

91.534 

-1.5384615 

126 

53.3 

84 

92.057 

-1.3 

2 

92.707 

-1.5384615 

128 

54 

82.7 

89.584 

-1.3 

2 

90.234 

-1.5384615 

128 

54.1 

82.8 

90.857 

-1.2 

2 

91.457 

-1.6666667 

130 

54.8 

81.5 

88.384 

-1.2 

2 

88.984 

-1.6666667 

130 

54.7 

81.4 

89.457 

-1.4 

2 

90.157 

-1.4285714 

132 

55.9 

80.6 

87.484 

-0.9 

2 

87.934 

-2.2222222 

132 

55.4 

80.1 

88.157 

-1.3 

2 

88.807 

-1.5384615 

134 

56.8 

79.5 

86.384 

-1.1 

2 

86.934 

-1.8181818 

134 

56.9 

79.6 

87.657 

-0.5 

2 

87.907 

-4 

136 

58 

78.7 

85.584 

-0.8 

2 

85.984 

-2.5 

136 

57.9 

78.6 

86.657 

-1 

2 

87.157 

-2 

138 

59.1 

77.8 

84.684 

-0.9 

2 

85.134 

-2.2222222 

138 

59.5 

78.2 

86.257 

-0.4 

2 

86.457 

-5 

140 

60.2 

76.9 

83.784 

-0.9 

2 

84.234 

-2.2222222 

140 

60.4 

77.1 

85.157 

-1.1 

2 

85.707 

-1.8181818 

144 

63.3 

76 

82.884 

-0.9 

4 

83.334 

-4.4444444 

144 

62.6 

75.3 

83.357 

-1.8 

4 

84.257 

-2.2222222 

148 

65.8 

74.5 

81.384 

-1.5 

4 

82.134 

-2.6666667 

148 

65.4 

74.1 

82.157 

-1.2 

4 

82.757 

-3.3333333 

Data  from  the  Camera  Objective  at  370  nm 


Forward 

Reverse 

Si 

Sq 

ASn 

ASi 

Si 

dnai 

S_Q 

ASn 

ASi 

80 

88.2 

164.9 

171.784 

n/a 

n/a 

n/a 

n/a 

80 

87.4 

164.1 

172.157 

n/a 

n/a 

n/a 

n/a 

81 

82.7 

158.4 

165.284 

-6.5 

1 

168.534 

-0.1538462 

81 

84 

159.7 

167.757 

A.4 

1 

169.957 

-0.2272727 

82 

78.1 

152.8 

159.684 

-5.6 

1 

162.484 

-0.1785714 

82 

80.8 

155.5 

163.557 

■4.2 

1 

165.657 

-0.2380952 

83 

75.4 

149.1 

155.984 

-3.7 

1 

157.834 

-0.2702703 

83 

77.6 

151.3 

159.357 

4.2 

1 

161.457 

-0.2380952 

84 

72.9 

145.6 

152.484 

-3.5 

1 

154.234 

-0.2857143 

84 

71.7 

144.4 

152.457 

-6.9 

1 

155.907 

-0.1449275 

85 

69.5 

141.2 

148.084 

-4.4 

1 

150.284 

-0.2272727 

85 

71.5 

143.2 

151.257 

-1.2 

1 

151.857 

-0.8333333 

86 

66.4 

137.1 

143.984 

-4.1 

1 

146.034 

-0.2439024 

86 

68.8 

139.5 

147.557 

-3.7 

1 

149.407 

-0.2702703 

87 

64.5 

134.2 

141.084 

-2.9 

1 

142.534 

-0.3448276 

87 

68.3 

138 

146.057 

-1.5 

1 

146.807 

-0.6666667 

88 

63.6 

132.3 

139.184 

-1.9 

1 

140.134 

-0.5263158 

88 

68 

136.7 

144.757 

-1.3 

1 

145.407 

-0.7692308 

89 

62.6 

130.3 

137.184 

-2 

1 

138.184 

-0.5 

89 

64.8 

132.5 

140.557 

•4.2 

1 

142.657 

-0.2380952 

90 

61.3 

128 

134.884 

-2.3 

1 

136.034 

-0.4347826 

90 

64 

130.7 

138.757 

-1.8 

1 

139.657 

-0.5555556 

91 

58.6 

124.3 

131.184 

-3.7 

1 

133.034 

-0.2702703 

91 

62.4 

128.1 

136.157 

-2.6 

1 

137.457 

-0.3846154 

92 

57.6 

122.3 

129.184 

-2 

1 

130.184 

-0.5 

92 

61.7 

126.4 

134.457 

-1.7 

1 

135.307 

-0.5882353 

93 

56.4 

120.1 

126.984 

-2.2 

1 

128.084 

-0.4545455 

93 

60.1 

123.8 

131.857 

-2.6 

1 

133.157 

-0.3846154 

94 

55.8 

118.5 

125.384 

-1.6 

1 

126.184 

-0.625 

94 

59.5 

122.2 

130.257 

-1.6 

1 

131.057 

-0.625 

95 

55.1 

116.8 

123.684 

-1.7 

1 

124.534 

-0.5882353 

95 

59.1 

120.8 

128.857 

-1.4 

1 

129.557 

-0.7142857 

96 

55 

115.7 

122.584 

-1.1 

1 

123.134 

-0.9090909 

96 

58.2 

118.9 

126.957 

-1.9 

1 

127.907 

-0.5263158 

97 

51.8 

111.5 

118.384 

-4.2 

1 

120.484 

-0.2380952 

97 

56.6 

116.3 

124.357 

-2.6 

1 

125.657 

-0.3846154 

98 

51.6 

110.3 

117.184 

-1.2 

1 

117.784 

-0.8333333 

98 

55.2 

113.9 

121.957 

-2.4 

1 

123.157 

-0.4166667 

99 

51.9 

109.6 

116.484 

-0.7 

1 

116.834 

-1.4285714 

99 

55.3 

113 

121.057 

-0.9 

1 

121.507 

-1.1111111 

100 

52.4 

109.1 

115.984 

-0.5 

1 

116.234 

-2 

100 

54 

110.7 

118.757 

-2.3 

1 

119.907 

-0.4347826 

102 

51.7 

106.4 

113.284 

-2.7 

2 

114.634 

-0.7407407 

CM 

O 

53.1 

107.8 

115.857 

-2.9 

2 

117.307 

-0.6896552 

104 

51.8 

104.5 

111.384 

-1.9 

2 

112.334 

-1.0526316 

104 

52.5 

105.2 

113.257 

-2.6 

2 

114.557 

-0.7692308 

106 

51.9 

102.6 

109.484 

-1.9 

2 

110.434 

-1.0526316 

106 

49.2 

99.9 

107.957 

-5.3 

2 

110.607 

-0.3773585 

108 

51.8 

100.5 

107.384 

-2.1 

2 

108.434 

-0.952381 

o 

CO 

51 

99.7 

107.757 

-0.2 

2 

107.857 

-10 

110 

52.4 

99.1 

105.984 

-1.4 

2 

106.684 

-1.4285714 

110 

52.4 

99.1 

107.157 

-0.6 

2 

107.457 

-3.3333333 

112 

53 

97.7 

104.584 

-1.4 

2 

105.284 

-1.4285714 

112 

52.2 

96.9 

104.957 

-2.2 

2 

106.057 

-0.9090909 

114 

51.9 

94.6 

101.484 

-3.1 

2 

103.034 

-0.6451613 

114 

53.2 

95.9 

103.957 

-1 

2 

104.457 

-2 

116 

52.6 

93.3 

100.184 

-1.3 

2 

100.834 

-1.5384615 

116 

54.2 

94.9 

102.957 

-1 

2 

103.457 

-2 

118 

53.2 

91.9 

98.784 

-1.4 

2 

99.484 

-1.4285714 

118 

55.8 

94.5 

102.557 

-0.4 

2 

102.757 

-5 

120 

54 

90.7 

97.584 

-1.2 

2 

98.184 

-1.6666667 

120 

57.6 

94.3 

102.357 

-0.2 

2 

102.457 

-10 

122 

54.6 

89.3 

96.184 

-1.4 

2 

96.884 

-1.4285714 

122 

59.5 

94.2 

102.257 

-0.1 

2 

102.307 

-20 

124 

55.5 

88.2 

95.084 

-1.1 

2 

95.634 

-1.8181818 

124 

61.4 

94.1 

102.157 

-0.1 

2 

102.207 

-20 

126 

56.4 

87.1 

93.984 

-1.1 

2 

94.534 

-1.8181818 

CO 

CM 

63 

93.7 

101.757 

-0.4 

2 

101.957 

-5 

128 

57.8 

86.5 

93.384 

-0.6 

2 

93.684 

-3.3333333 

CO 

CM 

63.8 

92.5 

100.557 

-1.2 

2 

101.157 

-1.6666667 

130 

58.9 

85.6 

92.484 

-0.9 

2 

92.934 

-2.2222222 

130 

64.4 

91.1 

99.157 

-1.4 

2 

99.857 

-1.4285714 

132 

60.2 

84.9 

91.784 

-0.7 

2 

92.134 

-2.8571429 

132 

64.7 

89.4 

97.457 

-1.7 

2 

98.307 

-1.1764706 

134 

61.6 

84.3 

91.184 

-0.6 

2 

91.484 

-3.3333333 

134 

66.4 

89.1 

97.157 

-0.3 

2 

97.307 

-6.6666667 

136 

63 

83.7 

90.584 

-0.6 

2 

90.884 

-3.3333333 

136 

67.2 

87.9 

95.957 

-1.2 

2 

96.557 

-1.6666667 

138 

64.5 

83.2 

90.084 

-0.5 

2 

90.334 

-4 

138 

68.4 

87.1 

95.157 

-0.8 

2 

95.557 

-2.5 

140 

65.4 

82.1 

88.984 

-1.1 

2 

89.534 

-1.8181818 

140 

70.3 

87 

95.057 

-0.1 

2 

95.107 

-20 

144 

67.3 

80 

86.884 

-2.1 

4 

87.934 

-1.9047619 

144 

72.1 

84.8 

92.857 

-2.2 

4 

93.957 

-1.8181818 

148 

69.9 

78.6 

85.484 

-1.4 

4 

86.184 

-2.8571429 

CO 

73 

81.7 

89.757 

-3.1 

4 

91.307 

-1.2903226 

f^FiltemdDatafixiria 


Rmad 


255 

14.65 

157 

rVa 

rVa 

rVa 

25 

14.8 

16 

-0.5 

■>151 

0.8 

24.5 

14.9 

16.1 

-0.5 

KO 

0.1 

24 

15.1 

1635 

-0.5 

0.25 

23.5 

154 

16.6 

-0.5 

mum 

0.25 

23 

15.55 

1685 

-0.5 

0.25 

225 

158 

17.1 

-0.5 

■igii 

0.25 

22 

16 

17.45 

-0.5 

0.2 

0.35 

21.5 

16.25 

17.6 

-0.5 

0.25 

0.15 

21 

16.56 

17.9 

-0.5 

0.8 

0.3 

23.5 

168 

1825 

-0.5 

0.25 

0.35 

23 

17.1 

1856 

-0.5 

0.8 

0.3 

19.5 

17.5 

19.05 

-0.5 

0.4 

0.5 

19 

17.8 

19.45 

-0.5 

0.8 

0.4 

185 

182 

19.95 

-0.5 

0.4 

0.5 

18 

188 

23.5 

-0.5 

0.6 

0.56 

17.5 

19.8 

21.1 

-0.5 

0.5 

0.6 

17 

19.8 

21.7 

-0.5 

0.5 

0.6 

165 

20.5 

224 

-0.5 

0.7 

0.7 

16 

21.2 

2845 

-0.5 

0.7 

1.05 

155 

22 

24.2 

-0.5 

0.8 

0.75 

15 

229 

258 

-0.5 

0.9 

1.1 

14.5 

289 

20.7 

-0.5 

1 

1.4 

14 

25.1 

27.9 

-0.5 

1.2 

1.2 

185 

^.75 

29.6 

-0.5 

1.65 

1.7 

18 

^.5 

81.5 

-0.5 

1.75 

1.9 

(/HI  dst3noesaeg\^incm) 


rVa 


A^90 

rVa 

2525 

24.75 

24.25 

23.75 

23.25 
2275 
2225 

21.75 

21.25 

23.75 

23.25 

19.75 

19.25 
1875 
1825 

17.75 

17.25 
1675 
1625 
1575 
1525 

14.75 

14.25 
1875 
1825 


rVa 

rVa 

26 

1875 

14.2 

rVa 

-0.3 

-0.6 

255 

18.9 

14.95 

-0.5 

-0.2 

-0.2 

25 

14.1 

152 

-0.5 

-0.4 

-0.5 

24.5 

14.8 

158 

-0.5 

-0.6 

-0.5 

24 

14.4 

155 

-0.5 

-0.3 

-0.5 

285 

14.7 

159 

-0.5 

-0.5 

-0.5 

23 

14.9 

162 

-0.5 

-0.4 

-0.7 

225 

15.2 

164 

-0.5 

-0.5 

-0.3 

22 

15.5 

168 

-0.5 

-0.6 

-0.6 

21.5 

15.8 

17.2 

-0.5 

-0.5 

-0.7 

21 

16.1 

17.5 

-0.5 

-0.6 

-0.6 

20.5 

165 

18 

-0.5 

-0.8 

-1 

20 

169 

1845 

-0.5 

-0.6 

-0.8 

19.5 

17.4 

19.05 

-0.5 

-0.8 

-1 

19 

17.8 

19.5 

-0.5 

-1.2 

-1.1 

185 

18.4 

20.2 

-0.5 

-1 

-1.2 

18 

19 

21 

-0.5 

-1 

-1.2 

17.5 

19.55 

21.8 

-0.5 

-1.4 

-1.4 

17 

23.2 

226 

-0.5 

-1.4 

-21 

165 

21.1 

288 

-0.5 

-1.6 

-1.5 

16 

22 

24.95 

-0.5 

-1.8 

-22 

155 

23l05 

268 

-0.5 

-2 

-28 

15 

24.4 

27.8 

-0.5 

-24 

-24 

14.5 

258 

S.66 

-0.5 

-83 

-3.4 

14 

27.9 

328 

-0.5 

-3.5 

-88 

185 

29.7 

35 

-0.5 

0.55 

0.65 

0.9 

0.9 

1.05 

1.35 

1.4 

21 

1.8 


Ad  90 

rVa 

rVa 

rVa 

rVa 

0.75 

2575 

-0.8 

-1.5 

0.25 

2525 

-0.4 

-0.5 

0.1 

24.75 

-0.4 

-0.2 

0.2 

24.25 

-0.2 

-0.4 

0.4 

2875 

-0.6 

-0.8 

0.3 

2825 

-0.4 

-0.6 

0.2 

2275 

-0.6 

-0.4 

0.4 

2225 

-0.6 

-0.8 

0.4 

21.75 

-0.6 

-0.8 

0.3 

21.25 

-0.6 

-0.6 

0.5 

20.75 

-0.8 

-1 

0.45 

20.25 

-0.8 

-0.9 

0.6 

19.75 

-1 

-1.2 

0.45 

19.25 

-0.8 

-0.9 

0.7 

1875 

-1.2 

-1.4 

0.8 

1825 

-1.2 

-1.6 

0.8 

17.75 

-1.1 

-1.6 

0.8 

17.25 

-1.8 

-1.6 

1.2 

1675 

-1.8 

-24 

1.15 

1625 

-1.8 

-28 

1.35 

1575 

-21 

-27 

1.5 

1525 

-27 

-8 

1.85 

14.75 

-28 

-87 

265 

14.25 

4.2 

-58 

27 

1875 

-8.6 

-54 

OnpacLsiB 


Rjmd 

FbBse 

l^s 

Aqsd 

23 

155 

168 

da 

da 

da 

da 

da 

da 

23 

15 

162 

da 

da 

da 

da 

da 

da 

225 

157 

17.1 

-Q5 

02 

03 

2275 

-04 

-06 

225 

1515 

164 

-05 

015 

02 

22® 

-03 

-04 

22 

16 

17.35 

-Q5 

03 

025 

2225 

-06 

-05 

22 

154 

168 

-05 

0® 

04 

22® 

-05 

-08 

21.5 

163 

17.7 

-Q5 

03 

0$ 

21.75 

-06 

-07 

21.5 

158 

172 

-05 

04 

04 

21.® 

-08 

-08 

21 

1655 

18 

-Q5 

025 

03 

2125 

-05 

-06 

21 

1615 

176 

-05 

0® 

04 

21.® 

-07 

-08 

235 

168 

183 

-Q5 

025 

03 

2D75 

-05 

-06 

2D5 

1645 

17.® 

-05 

03 

0® 

®® 

-06 

-07 

2D 

1715 

187 

-Q5 

0$ 

04 

2D25 

-07 

-08 

2D 

168 

184 

-05 

0® 

045 

®® 

-07 

-09 

195 

175 

19C5 

-Q5 

0$ 

0$ 

1975 

-07 

-07 

195 

17.3 

18® 

-05 

05 

0® 

19® 

-1 

-1.1 

19 

179 

1955 

-Q5 

04 

05 

1925 

-08 

-1 

19 

1775 

2D® 

-05 

045 

1.6 

19® 

-09 

-32 

185 

183 

19S5 

-Q5 

04 

04 

1875 

-08 

-08 

185 

183 

2D1 

-05 

0® 

-045 

18® 

-1.1 

09 

18 

188 

2D$ 

-Q5 

05 

06 

1825 

-1 

-1.2 

18 

188 

21 

-05 

05 

09 

18® 

-1 

-1.8 

175 

192 

21.15 

-Q5 

04 

06 

17.75 

-08 

-12 

17.5 

195 

21.7 

-05 

07 

07 

17® 

-1.4 

-1.4 

17 

1985 

21.7 

-Q5 

0® 

0$ 

1725 

-1.3 

-1.1 

17 

2D1 

226 

-05 

06 

09 

17® 

-1.2 

-1.8 

165 

2D5 

2245 

-Q5 

0$ 

075 

1675 

-1.3 

-1.5 

165 

21.® 

236 

-05 

0® 

1 

16® 

-1.9 

-2 

16 

21.3 

232 

-Q5 

08 

075 

1625 

-1.6 
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APPENDIX  B.  MATLAB  CODE 


%Calculatioii  of  matrix  elements  for  the 
%Primary  Objective  from  data  collected  using 
%a  10  nm  bandpass  filter  at  220  nm 

clear  all 

%This  portion  of  the  program  will  load 
%the  image  and  object  distance  measurements 
%taken  in  the  lab,  fit  them  to  a  general 
%hyperbolic  function,  and  use  that  function 
%to  generate  another  set  of  data  for  further 
%analysis 

%Import  the  initial  data 
load  Lens3_220nmdata.dat; 

%Create  vectors  for  analysis  with  the  curve 
%fitting  function 

fwdxdata  =[[1  0  0  0]*Lens3_220nmdata']; 
fwdydata  =[[0  1  0  0]*Lens3_220nmdata']; 
revxdata  =[[0  0  1  0]*Lens3_220nmdata']; 
revydata  =[[0  0  0  l]*Lens3_220nmdata']; 

options=optimset(']VtaxFunEvals',le30,'maxiter',5000); 

%Use  f  s  for  initial  curve  fit  of  forward  data 
fD=[l  1  1  1]; 

[fresnorm]=lsqcurvefit(@GenHyperbola,fO, fwdxdata,  fwdydata); 

%Use  r's  for  initial  curve  fit  of  reverse  data 
r0=[l  1  1  1]; 

[r  resnorm]=lsqcurvefit(@GenHyperbola,rO,revxdata,revydata); 

%Generation  of  new  image  and  object  distance  data 
fwdsi=35:0.1:100; 

fwdso=((f(  1 ) .  *fwdsi+f(2))  ./(f(3) .  *  fwdsi+f(4))) ; 
revsi=35:0. 1:100; 

revso=((r(  1 ).  *revsi+r(2))./(r(3 ) .  *revsi+r(4))) ; 

%Plot  results  for  this  section 

subp  lot(2 ,2 , 1 ) ;  plot  (fwdxdata,  fwdydata,  'ro ') ; 
title('Forward  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
hold; 

plot(fwdsi,fwdso,'b'); 
ylim([0  150]); 

subp  lot(2,2,3);p  lot(revxdata,re  vydata,  'ro ') ; 
title('Reverse  Data'); 
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xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
hold; 

plot(revsi,revso,'b'); 
ylim([0  150]); 

%Generation  of  new  dsi/dso  data 

fwddsidso=(f(l)-*f(4)-f(2).*f(3))./(f(3).*fwdso-f(l)).^2; 

revdsidso=(r(l).*r(4)-r(2).*r(3))./(r(3).*revso-r(l)).^2; 

%This  portion  of  the  program  will  utilize  the  data 
%generated  earlier  to  calculate  the  matrix  elements 
%of  the  complex  lens  system 

%Use  a's  for  forward  data 
a0=[l  1]; 

[a  resnorm]=lsqcurvefit(@thindsidso2,a0,fwdso,fwddsidso) 

%Use  b's  for  reverse  data 
b0=[l  1]; 

[b  resnorm]=lsqcurvefit(@thindsidso2,b0,revso,revdsidso) 

%Calculate  Lens  Matrix  Elements 

mc=-l/((b(l)+a(l)y2); 

md=mc*(-((a(l)+b(l))/2)-a(2)); 

ma=l-mc*b(2); 

mb=(ma  *  md- 1 )/ me ; 

lensmatrix=[ma  mb;mc  md] 

mp=md/mc; 

mq=-ma/mc; 

fwdfocallength=abs(mq) 

revfocallength=abs(mp) 

xdatum=0:0.1:90; 

fwdydatum=(-(a(  1 )  .^2)  ./(xdatum-a(2)-a(  1 ))  .^2) ; 
revydatum=(-(b(l).^2)./(xdatum-b(2)-b(l)).^2); 
subplot(2 , 2, 2) ;  p lot(xdatum, fwdydatum) ; 
ylim([-50,0]); 

title('Forward  Data  (Asymptote  is  reverse  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

plot(fwdso,fwddsidso,'ro'); 

xlim([0,100]); 

subplot(2 , 2, 4) ;  p lot(xdatum,revy datum) ; 
ylim([-50,0]); 

title('Reverse  Data  (Asymptote  is  forward  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

p  lot(re  vso  ,revdsidso ,  'ro ') ; 
xlim([0,100]); 
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%Calculatioii  of  matrix  elements  for  the 
%Primary  Objective  from  data  collected  using 
%a  10  nm  bandpass  filter  at  300  nm 

clear  all 

%This  portion  of  the  program  will  load 
%the  image  and  object  distance  measurements 
%taken  in  the  lab,  fit  them  to  a  general 
%hyperbolic  function,  and  use  that  function 
%to  generate  another  set  of  data  for  further 
%analysis 

%Import  the  initial  data 

load  Lens3_300nmdata.dat; 

%Create  vectors  for  analysis  with  the  curve 
%fitting  function 

fwdxdata  =[[1  0  0  0]*Lens3_300nmdata']; 
fwdydata  =[[0  1  0  0]*Lens3_300nmdata']; 
revxdata  =[[0  0  1  0]*Lens3_300nmdata']; 
revydata  =[[0  0  0  l]*Lens3_300nmdata']; 

options=optimset(']VtaxFunEvals',le30,'maxiter',5000); 

%Use  f  s  for  initial  curve  fit  of  forward  data 
fD=[l  1  1  1]; 

[fresnorm]=lsqcurvefit(@GenHyperbola,fO, fwdxdata,  fwdydata); 

%Use  r's  for  initial  curve  fit  of  reverse  data 
r0=[l  1  1  1]; 

[r  resnorm]=lsqcurvefit(@GenHyperbola,rO,revxdata,revydata); 
%Generation  of  new  image  and  object  distance  data 
fwdsi=35:0.1:100; 

fwdso=((f(  1 ) .  *fwdsi+f(2))  ./(f(3) .  *  fwdsi+f(4))) ; 
revsi=35:0. 1:100; 

revso=((r(  1 ).  *revsi+r(2))./(r(3 ) .  *revsi+r(4))) ; 

%Plot  results  for  this  section 

subp  lot(2,2,l);p  lot  (fwdxdata,  fwdydata,  'ro ') ; 
title('Forward  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
hold; 

plot(fwdsi,fwdso,'b'); 
ylim([0  150]); 

subplot(2 ,2 , 3 ) ;  plot(revxdata,re vy data, 'ro ') ; 
title('Reverse  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
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hold; 

plot(revsi,revso,'b'); 
ylim([0  150]); 

%Generation  of  new  dsi/dso  data 

fwddsidso=(f(l)-*f(4)-f(2).*f(3))./(f(3).*fwdso-f(l)).^2; 

revdsidso=(r(l).*r(4)-r(2).*r(3))./(r(3).*revso-r(l)).^2; 

%This  portion  of  the  program  will  utilize  the  data 
%generated  earlier  to  caleulate  the  matrix  elements 
%of  the  complex  lens  system 

%Use  a's  for  forward  data 
a0=[l  1]; 

[a  resnorm]=lsqcurvefit(@thindsidso2,a0,fwdso,fwddsidso) 

%Use  b's  for  reverse  data 
b0=[l  1]; 

[b  resnorm]=lsqcurvefit(@thindsidso2,b0,revso,revdsidso) 

%Calculate  Lens  Matrix  Elements 

mc=-l/((b(l)+a(l)y2); 

md=mc*(-((a(l)+b(l))/2)-a(2)); 

ma=l-mc*b(2); 

mb=(ma  *  md- 1 )/ me ; 

lensmatrix=[ma  mb;mc  md] 

mp=md/mc; 

mq=-ma/mc; 

fwdfocallength=abs(mq) 

revfocallength=abs(mp) 

xdatum=0:0.1:90; 

fwdydatum=(-(a(  1 )  .^2)  ./(xdatum-a(2)-a(  1 ))  .^2) ; 
revydatum=(-(b(l).^2)./(xdatum-b(2)-b(l)).^2); 
subplot(2,2,2);  plot(xdatum,fwdydatum); 
ylim([-50,0]); 

title('Forward  Data  (Asymptote  is  reverse  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

plot(fwdso,fwddsidso,'ro'); 

xlim([0,100]); 

subplot(2 , 2, 4) ;  p lot(xdatum,revy datum) ; 
ylim([-50,0]); 

title('Reverse  Data  (Asymptote  is  forward  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

p  lot(re  vso  ,revdsidso ,  'ro ') ; 
xlim([0,100]); 
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%Calculatioii  of  matrix  elements  for  the 
%Primary  Objective  from  data  collected  using 
%a  10  nm  bandpass  filter  at  334  nm 

clear  all 

%This  portion  of  the  program  will  load 
%the  image  and  object  distance  measurements 
%taken  in  the  lab,  fit  them  to  a  general 
%hyperbolic  function,  and  use  that  function 
%to  generate  another  set  of  data  for  further 
%analysis 

%Import  the  initial  data 

load  Lens3_334nmdata.dat; 

%Create  vectors  for  analysis  with  the  curve 
%fitting  function 

fwdxdata  =[[1  0  0  0]*Lens3_334nmdata']; 
fwdydata  =[[0  1  0  0]*Lens3_334nmdata']; 
revxdata  =[[0  0  1  0]*Lens3_334nmdata']; 
revydata  =[[0  0  0  l]*Lens3_334nmdata']; 

options=optimset(']VtaxFunEvals',le30,'maxiter',5000); 

%Use  f  s  for  initial  curve  fit  of  forward  data 
fD=[l  1  1  1]; 

[fresnorm]=lsqcurvefit(@GenHyperbola,fO, fwdxdata,  fwdydata); 

%Use  r's  for  initial  curve  fit  of  reverse  data 
r0=[l  1  1  1]; 

[r  resnorm]=lsqcurvefit(@GenHyperbola,rO,revxdata,revydata); 
%Generation  of  new  image  and  object  distance  data 
fwdsi=35:0.1:100; 

fwdso=((f(  1 ) .  *fwdsi+f(2))  ./(f(3) .  *  fwdsi+f(4))) ; 
revsi=35:0. 1:100; 

revso=((r(  1 ).  *revsi+r(2))./(r(3 ) .  *revsi+r(4))) ; 

%Plot  results  for  this  section 

subp  lot(2,2,l);p  lot  (fwdxdata,  fwdydata,  'ro ') ; 
title('Forward  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
hold; 

plot(fwdsi,fwdso,'b'); 
ylim([0  150]); 

subplot(2 ,2 , 3 ) ;  plot(revxdata,re vy data, 'ro ') ; 
title('Reverse  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
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hold; 

plot(revsi,revso,'b'); 
ylim([0  150]); 

%Generation  of  new  dsi/dso  data 

fwddsidso=(f(l)-*f(4)-f(2).*f(3))./(f(3).*fwdso-f(l)).^2; 

revdsidso=(r(l).*r(4)-r(2).*r(3))./(r(3).*revso-r(l)).^2; 

%This  portion  of  the  program  will  utilize  the  data 
%generated  earlier  to  caleulate  the  matrix  elements 
%of  the  complex  lens  system 

%Use  a's  for  forward  data 
a0=[l  1]; 

[a  resnorm]=lsqcurvefit(@thindsidso2,a0,fwdso,fwddsidso) 

%Use  b's  for  reverse  data 
b0=[l  1]; 

[b  resnorm]=lsqcurvefit(@thindsidso2,b0,revso,revdsidso) 

%Calculate  Lens  Matrix  Elements 

mc=-l/((b(l)+a(l)y2); 

md=mc*(-((a(l)+b(l))/2)-a(2)); 

ma=l-mc*b(2); 

mb=(ma  *  md- 1 )/ me ; 

lensmatrix=[ma  mb;mc  md] 

mp=md/mc; 

mq=-ma/mc; 

fwdfocallength=abs(mq) 

revfocallength=abs(mp) 

xdatum=0:0.1:90; 

fwdydatum=(-(a(  1 )  .^2)  ./(xdatum-a(2)-a(  1 ))  .^2) ; 
revydatum=(-(b(l).^2)./(xdatum-b(2)-b(l)).^2); 
subplot(2,2,2);  plot(xdatum,fwdydatum); 
ylim([-50,0]); 

title('Forward  Data  (Asymptote  is  reverse  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

plot(fwdso,fwddsidso,'ro'); 

xlim([0,100]); 

subplot(2 , 2, 4) ;  p lot(xdatum,revy datum) ; 
ylim([-50,0]); 

title('Reverse  Data  (Asymptote  is  forward  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

p  lot(re  vso  ,revdsidso ,  'ro ') ; 
xlim([0,100]); 
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%Calculatioii  of  matrix  elements  for  the 
%Primary  Objective  from  data  collected  using 
%a  10  nm  bandpass  filter  at  370  nm 

clear  all 

%This  portion  of  the  program  will  load 
%the  image  and  object  distance  measurements 
%taken  in  the  lab,  fit  them  to  a  general 
%hyperbolic  function,  and  use  that  function 
%to  generate  another  set  of  data  for  further 
%analysis 

%Import  the  initial  data 

load  Lens3_370nmdata.dat; 

%Create  vectors  for  analysis  with  the  curve 
%fitting  function 

fwdxdata  =[[1  0  0  0]*Lens3_370nmdata']; 
fwdydata  =[[0  1  0  0]*Lens3_370nmdata']; 
revxdata  =[[0  0  1  0]*Lens3_370nmdata']; 
revydata  =[[0  0  0  l]*Lens3_370nmdata']; 

options=optimset(']VtaxFunEvals',le30,'maxiter',5000); 

%Use  f  s  for  initial  curve  fit  of  forward  data 
fD=[l  1  1  1]; 

[fresnorm]=lsqcurvefit(@GenHyperbola,fO, fwdxdata,  fwdydata); 

%Use  r's  for  initial  curve  fit  of  reverse  data 
r0=[l  1  1  1]; 

[r  resnorm]=lsqcurvefit(@GenHyperbola,rO,revxdata,revydata); 
%Generation  of  new  image  and  object  distance  data 
fwdsi=35:0.1:100; 

fwdso=((f(  1 ) .  *fwdsi+f(2))  ./(f(3) .  *  fwdsi+f(4))) ; 
revsi=35:0. 1:100; 

revso=((r(  1 ).  *revsi+r(2))./(r(3 ) .  *revsi+r(4))) ; 

%Plot  results  for  this  section 

subp  lot(2,2,l);p  lot  (fwdxdata,  fwdydata,  'ro ') ; 
title('Forward  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
hold; 

plot(fwdsi,fwdso,'b'); 
ylim([0  150]); 

subplot(2 ,2 , 3 ) ;  plot(revxdata,re vy data, 'ro ') ; 
title('Reverse  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
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hold; 

plot(revsi,revso,'b'); 
ylim([0  150]); 

%Generation  of  new  dsi/dso  data 

fwddsidso=(f(l)-*f(4)-f(2).*f(3))./(f(3).*fwdso-f(l)).^2; 

revdsidso=(r(l).*r(4)-r(2).*r(3))./(r(3).*revso-r(l)).^2; 

%This  portion  of  the  program  will  utilize  the  data 
%generated  earlier  to  caleulate  the  matrix  elements 
%of  the  complex  lens  system 

%Use  a's  for  forward  data 
a0=[l  1]; 

[a  resnorm]=lsqcurvefit(@thindsidso2,a0,fwdso,fwddsidso) 

%Use  b's  for  reverse  data 
b0=[l  1]; 

[b  resnorm]=lsqcurvefit(@thindsidso2,b0,revso,revdsidso) 

%Calculate  Lens  Matrix  Elements 

mc=-l/((b(l)+a(l)y2); 

md=mc*(-((a(l)+b(l))/2)-a(2)); 

ma=l-mc*b(2); 

mb=(ma  *  md- 1 )/ me ; 

lensmatrix=[ma  mb;mc  md] 

mp=md/mc; 

mq=-ma/mc; 

fwdfocallength=abs(mq) 

revfocallength=abs(mp) 

xdatum=0:0.1:90; 

fwdydatum=(-(a(  1 )  .^2)  ./(xdatum-a(2)-a(  1 ))  .^2) ; 
revydatum=(-(b(l).^2)./(xdatum-b(2)-b(l)).^2); 
subplot(2,2,2);  plot(xdatum,fwdydatum); 
ylim([-50,0]); 

title('Forward  Data  (Asymptote  is  reverse  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

plot(fwdso,fwddsidso,'ro'); 

xlim([0,100]); 

subplot(2 , 2, 4) ;  p lot(xdatum,revy datum) ; 
ylim([-50,0]); 

title('Reverse  Data  (Asymptote  is  forward  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

p  lot(re  vso  ,revdsidso ,  'ro ') ; 
xlim([0,100]); 
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%Calculatioii  of  matrix  elements  for  the 
%Camera  Objective  from  data  collected  using 
%a  10  nm  bandpass  filter  at  220  nm 

clear  all 

%This  portion  of  the  program  will  load 
%the  image  and  object  distance  measurements 
%taken  in  the  lab,  fit  them  to  a  general 
%hyperbolic  function,  and  use  that  function 
%to  generate  another  set  of  data  for  further 
%analysis 

%Import  the  initial  data 

load  Lensl_220nmdata.dat; 

%Create  vectors  for  analysis  with  the  curve 
%fitting  function 

fwdxdata  =[[1  0  0  0]*Lensl_220nmdata']; 
fwdydata  =[[0  1  0  0]*Lensl_220nmdata']; 
revxdata  =[[0  0  1  0]*Lensl_220nmdata']; 
revydata  =[[0  0  0  l]*Lensl_220nmdata']; 

options=optimset(']VtaxFunEvals',le30,'maxiter',5000); 

%Use  f  s  for  initial  curve  fit  of  forward  data 
fD=[l  1  1  1]; 

[fresnorm]=lsqcurvefit(@GenHyperbola,fO, fwdxdata,  fwdydata); 

%Use  r's  for  initial  curve  fit  of  reverse  data 
r0=[l  1  1  1]; 

[r  resnorm]=lsqcurvefit(@GenHyperbola,rO,revxdata,revydata); 
%Generation  of  new  image  and  object  distance  data 
fwdsi=60:0. 1:150; 

fwdso=((f(  1 ) .  *fwdsi+f(2))  ./(f(3) .  *  fwdsi+f(4))) ; 
revsi=60:0. 1:150; 

revso=((r(  1 ).  *revsi+r(2))./(r(3 ) .  *revsi+r(4))) ; 

%Plot  results  for  this  section 

subp  lot(2,2,l);p  lot  (fwdxdata,  fwdydata,  'ro ') ; 
title('Forward  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
hold; 

plot(fwdsi,fwdso,'b'); 
ylim([0  200]); 

subplot(2 ,2 , 3 ) ;  plot(revxdata,re vy data, 'ro ') ; 
title('Reverse  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
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hold; 

plot(revsi,revso,'b'); 
ylim([0  200]); 

%Generation  of  new  dsi/dso  data 

fwddsidso=(f(l)-*f(4)-f(2).*f(3))./(f(3).*fwdso-f(l)).^2; 

revdsidso=(r(l).*r(4)-r(2).*r(3))./(r(3).*revso-r(l)).^2; 

%This  portion  of  the  program  will  utilize  the  data 
%generated  earlier  to  caleulate  the  matrix  elements 
%of  the  complex  lens  system 

%Use  a's  for  forward  data 
a0=[l  1]; 

[a  resnorm]=lsqcurvefit(@thindsidso2,a0,fwdso,fwddsidso) 

%Use  b's  for  reverse  data 
b0=[l  1]; 

[b  resnorm]=lsqcurvefit(@thindsidso2,b0,revso,revdsidso) 

%Calculate  Lens  Matrix  Elements 

mc=-l/((b(l)+a(l)y2); 

md=mc*(-((a(l)+b(l))/2)-a(2)); 

ma=l-mc*b(2); 

mb=(ma  *  md- 1 )/ me ; 

lensmatrix=[ma  mb;mc  md] 

mp=md/mc; 

mq=-ma/mc; 

fwdfocallength=abs(mq) 

revfocallength=abs(mp) 

xdatum=0:0.1:90; 

fwdydatum=(-(a(  1 )  .^2)  ./(xdatum-a(2)-a(  1 ))  .^2) ; 
revydatum=(-(b(l).^2)./(xdatum-b(2)-b(l)).^2); 
subplot(2,2,2);  plot(xdatum,fwdydatum); 
ylim([-50,0]); 

title('Forward  Data  (Asymptote  is  reverse  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

plot(fwdso,fwddsidso,'ro'); 

xlim([0,150]); 

subplot(2 , 2, 4) ;  p lot(xdatum,revy datum) ; 
ylim([-50,0]); 

title('Reverse  Data  (Asymptote  is  forward  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

p  lot(re  vso  ,revdsidso ,  'ro ') ; 
xlim([0,150]); 
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%Calculatioii  of  matrix  elements  for  the 
%Camera  Objective  from  data  collected  using 
%a  10  nm  bandpass  filter  at  300  nm 

clear  all 

%This  portion  of  the  program  will  load 
%the  image  and  object  distance  measurements 
%taken  in  the  lab,  fit  them  to  a  general 
%hyperbolic  function,  and  use  that  function 
%to  generate  another  set  of  data  for  further 
%analysis 

%Import  the  initial  data 

load  Lensl_300nmdata.dat; 

%Create  vectors  for  analysis  with  the  curve 
%fitting  function 

fwdxdata  =[[1  0  0  0]*Lensl_300nmdata']; 
fwdydata  =[[0  1  0  0]*Lensl_300nmdata']; 
revxdata  =[[0  0  1  0]*Lensl_300nmdata']; 
revydata  =[[0  0  0  l]*Lensl_300nmdata']; 

options=optimset(']VtaxFunEvals',le30,'maxiter',5000); 

%Use  f  s  for  initial  curve  fit  of  forward  data 
fD=[l  1  1  1]; 

[fresnorm]=lsqcurvefit(@GenHyperbola,fO, fwdxdata,  fwdydata); 

%Use  r's  for  initial  curve  fit  of  reverse  data 
r0=[l  1  1  1]; 

[r  resnorm]=lsqcurvefit(@GenHyperbola,rO,revxdata,revydata); 
%Generation  of  new  image  and  object  distance  data 
fwdsi=60:0. 1:150; 

fwdso=((f(  1 ) .  *fwdsi+f(2))  ./(f(3) .  *  fwdsi+f(4))) ; 
revsi=60:0. 1:150; 

revso=((r(  1 ).  *revsi+r(2))./(r(3 ) .  *revsi+r(4))) ; 

%Plot  results  for  this  section 

subp  lot(2,2,l);p  lot  (fwdxdata,  fwdydata,  'ro ') ; 
title('Forward  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
hold; 

plot(fwdsi,fwdso,'b'); 
ylim([0  200]); 

subplot(2 ,2 , 3 ) ;  plot(revxdata,re vy data, 'ro ') ; 
title('Reverse  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
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hold; 

plot(revsi,revso,'b'); 
ylim([0  200]); 

%Generation  of  new  dsi/dso  data 

fwddsidso=(f(l)-*f(4)-f(2).*f(3))./(f(3).*fwdso-f(l)).^2; 

revdsidso=(r(l).*r(4)-r(2).*r(3))./(r(3).*revso-r(l)).^2; 

%This  portion  of  the  program  will  utilize  the  data 
%generated  earlier  to  caleulate  the  matrix  elements 
%of  the  complex  lens  system 

%Use  a's  for  forward  data 
a0=[l  1]; 

[a  resnorm]=lsqcurvefit(@thindsidso2,a0,fwdso,fwddsidso) 

%Use  b's  for  reverse  data 
b0=[l  1]; 

[b  resnorm]=lsqcurvefit(@thindsidso2,b0,revso,revdsidso) 

%Calculate  Lens  Matrix  Elements 

mc=-l/((b(l)+a(l)y2); 

md=mc*(-((a(l)+b(l))/2)-a(2)); 

ma=l-mc*b(2); 

mb=(ma  *  md- 1 )/ me ; 

lensmatrix=[ma  mb;mc  md] 

mp=md/mc; 

mq=-ma/mc; 

fwdfocallength=abs(mq) 

revfocallength=abs(mp) 

xdatum=0:0.1:90; 

fwdydatum=(-(a(  1 )  .^2)  ./(xdatum-a(2)-a(  1 ))  .^2) ; 
revydatum=(-(b(l).^2)./(xdatum-b(2)-b(l)).^2); 
subplot(2,2,2);  plot(xdatum,fwdydatum); 
ylim([-50,0]); 

title('Forward  Data  (Asymptote  is  reverse  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

plot(fwdso,fwddsidso,'ro'); 

xlim([0,150]); 

subplot(2 , 2, 4) ;  p lot(xdatum,revy datum) ; 
ylim([-50,0]); 

title('Reverse  Data  (Asymptote  is  forward  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

p  lot(re  vso  ,revdsidso ,  'ro ') ; 
xlim([0,150]); 
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%Calculatioii  of  matrix  elements  for  the 
%Camera  Objective  from  data  collected  using 
%a  10  nm  bandpass  filter  at  334  nm 

clear  all 

%This  portion  of  the  program  will  load 
%the  image  and  object  distance  measurements 
%taken  in  the  lab,  fit  them  to  a  general 
%hyperbolic  function,  and  use  that  function 
%to  generate  another  set  of  data  for  further 
%analysis 

%Import  the  initial  data 

load  Lensl_334nmdata.dat; 

%Create  vectors  for  analysis  with  the  curve 
%fitting  function 

fwdxdata  =[[1  0  0  0]*Lensl_334nmdata']; 
fwdydata  =[[0  1  0  0]*Lensl_334nmdata']; 
revxdata  =[[0  0  1  0]*Lensl_334nmdata']; 
revydata  =[[0  0  0  l]*Lensl_334nmdata']; 

options=optimset(']VtaxFunEvals',le30,'maxiter',5000); 

%Use  f  s  for  initial  curve  fit  of  forward  data 
fD=[l  1  1  1]; 

[fresnorm]=lsqcurvefit(@GenHyperbola,fO, fwdxdata,  fwdydata); 

%Use  r's  for  initial  curve  fit  of  reverse  data 
r0=[l  1  1  1]; 

[r  resnorm]=lsqcurvefit(@GenHyperbola,rO,revxdata,revydata); 
%Generation  of  new  image  and  object  distance  data 
fwdsi=60:0. 1:150; 

fwdso=((f(  1 ) .  *fwdsi+f(2))  ./(f(3) .  *  fwdsi+f(4))) ; 
revsi=60:0. 1:150; 

revso=((r(  1 ).  *revsi+r(2))./(r(3 ) .  *revsi+r(4))) ; 

%Plot  results  for  this  section 

subp  lot(2,2,l);p  lot  (fwdxdata,  fwdydata,  'ro ') ; 
title('Forward  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
hold; 

plot(fwdsi,fwdso,'b'); 
ylim([0  200]); 

subplot(2 ,2 , 3 ) ;  plot(revxdata,re vy data, 'ro ') ; 
title('Reverse  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
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hold; 

plot(revsi,revso,'b'); 
ylim([0  200]); 

%Generation  of  new  dsi/dso  data 

fwddsidso=(f(l)-*f(4)-f(2).*f(3))./(f(3).*fwdso-f(l)).^2; 

revdsidso=(r(l).*r(4)-r(2).*r(3))./(r(3).*revso-r(l)).^2; 

%This  portion  of  the  program  will  utilize  the  data 
%generated  earlier  to  caleulate  the  matrix  elements 
%of  the  complex  lens  system 

%Use  a's  for  forward  data 
a0=[l  1]; 

[a  resnorm]=lsqcurvefit(@thindsidso2,a0,fwdso,fwddsidso) 

%Use  b's  for  reverse  data 
b0=[l  1]; 

[b  resnorm]=lsqcurvefit(@thindsidso2,b0,revso,revdsidso) 

%Calculate  Lens  Matrix  Elements 

mc=-l/((b(l)+a(l)y2); 

md=mc*(-((a(l)+b(l))/2)-a(2)); 

ma=l-mc*b(2); 

mb=(ma  *  md- 1 )/ me ; 

lensmatrix=[ma  mb;mc  md] 

mp=md/mc; 

mq=-ma/mc; 

fwdfocallength=abs(mq) 

revfocallength=abs(mp) 

xdatum=0:0.1:90; 

fwdydatum=(-(a(  1 )  .^2)  ./(xdatum-a(2)-a(  1 ))  .^2) ; 
revydatum=(-(b(l).^2)./(xdatum-b(2)-b(l)).^2); 
subplot(2,2,2);  plot(xdatum,fwdydatum); 
ylim([-50,0]); 

title('Forward  Data  (Asymptote  is  reverse  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

plot(fwdso,fwddsidso,'ro'); 

xlim([0,150]); 

subplot(2 , 2, 4) ;  p lot(xdatum,revy datum) ; 
ylim([-50,0]); 

title('Reverse  Data  (Asymptote  is  forward  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

p  lot(re  vso  ,revdsidso ,  'ro ') ; 
xlim([0,150]); 
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%Calculatioii  of  matrix  elements  for  the 
%Camera  Objective  from  data  collected  using 
%a  10  nm  bandpass  filter  at  370  nm 

clear  all 

%This  portion  of  the  program  will  load 
%the  image  and  object  distance  measurements 
%taken  in  the  lab,  fit  them  to  a  general 
%hyperbolic  function,  and  use  that  function 
%to  generate  another  set  of  data  for  further 
%analysis 

%Import  the  initial  data 

load  Lensl_370nmdata.dat; 

%Create  vectors  for  analysis  with  the  curve 
%fitting  function 

fwdxdata  =[[1  0  0  0]*Lensl_370nmdata']; 
fwdydata  =[[0  1  0  0]*Lensl_370nmdata']; 
revxdata  =[[0  0  1  0]*Lensl_370nmdata']; 
revydata  =[[0  0  0  l]*Lensl_370nmdata']; 

options=optimset(']VtaxFunEvals',le30,'maxiter',5000); 

%Use  f  s  for  initial  curve  fit  of  forward  data 
fD=[l  1  1  1]; 

[fresnorm]=lsqcurvefit(@GenHyperbola,fO, fwdxdata,  fwdydata); 

%Use  r's  for  initial  curve  fit  of  reverse  data 
r0=[l  1  1  1]; 

[r  resnorm]=lsqcurvefit(@GenHyperbola,rO,revxdata,revydata); 
%Generation  of  new  image  and  object  distance  data 
fwdsi=60:0. 1:150; 

fwdso=((f(  1 ) .  *fwdsi+f(2))  ./(f(3) .  *  fwdsi+f(4))) ; 
revsi=60:0. 1:150; 

revso=((r(  1 ).  *revsi+r(2))./(r(3 ) .  *revsi+r(4))) ; 

%Plot  results  for  this  section 

subp  lot(2,2,l);p  lot  (fwdxdata,  fwdydata,  'ro ') ; 
title('Forward  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
hold; 

plot(fwdsi,fwdso,'b'); 
ylim([0  200]); 

subplot(2 ,2 , 3 ) ;  plot(revxdata,re vy data, 'ro ') ; 
title('Reverse  Data'); 
xlabel('Object  Distance  (cm)'); 
ylabel('Image  Distance  (cm)'); 
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hold; 

plot(revsi,revso,'b'); 
ylim([0  200]); 

%Generation  of  new  dsi/dso  data 

fwddsidso=(f(l)-*f(4)-f(2).*f(3))./(f(3).*fwdso-f(l)).^2; 

revdsidso=(r(l).*r(4)-r(2).*r(3))./(r(3).*revso-r(l)).^2; 

%This  portion  of  the  program  will  utilize  the  data 
%generated  earlier  to  caleulate  the  matrix  elements 
%of  the  complex  lens  system 

%Use  a's  for  forward  data 
a0=[l  1]; 

[a  resnorm]=lsqcurvefit(@thindsidso2,a0,fwdso,fwddsidso) 

%Use  b's  for  reverse  data 
b0=[l  1]; 

[b  resnorm]=lsqcurvefit(@thindsidso2,b0,revso,revdsidso) 

%Calculate  Lens  Matrix  Elements 

mc=-l/((b(l)+a(l)y2); 

md=mc*(-((a(l)+b(l))/2)-a(2)); 

ma=l-mc*b(2); 

mb=(ma  *  md- 1 )/ me ; 

lensmatrix=[ma  mb;mc  md] 

mp=md/mc; 

mq=-ma/mc; 

fwdfocallength=abs(mq) 

revfocallength=abs(mp) 

xdatum=0:0.1:90; 

fwdydatum=(-(a(  1 )  .^2)  ./(xdatum-a(2)-a(  1 ))  .^2) ; 
revydatum=(-(b(l).^2)./(xdatum-b(2)-b(l)).^2); 
subplot(2,2,2);  plot(xdatum,fwdydatum); 
ylim([-50,0]); 

title('Forward  Data  (Asymptote  is  reverse  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

plot(fwdso,fwddsidso,'ro'); 

xlim([0,150]); 

subplot(2 , 2, 4) ;  p lot(xdatum,revy datum) ; 
ylim([-50,0]); 

title('Reverse  Data  (Asymptote  is  forward  focal  length)'); 

xlabel('Object  Distance  (cm)'); 

ylabel('dsi/dso'); 

hold; 

p  lot(re  vso  ,revdsidso ,  'ro ') ; 
xlim([0,150]); 
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%General  Hyperbolic  Function 

function  f=GenHyperbola(r,xdata) 
f^((r(  1 ) .  *xdata+r(2))./(r(3 ) .  *xdata+r(4))) ; 

%Function  relating  d$i/d$o  and  focal  length 
%for  a  theoretical  thin  lens 

function  f=thindsidso2(x,xdata) 

f^(-(x(  1 ) .  ^2)  ./(xdata-x(2)-x(  1 )) .  ^2) ; 
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